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Abstract 

Multi-particle form factors of local operators in integrable models in two dimensions 
seem to have the property that they factorize when one subset of the particles in 
the external states are boosted by a large rapidity with respect to the others. This 
remarkable property, which goes under the name of form factor clustering, was first 
observed by Smirnov in the 0(3) non-linear <r-model and has subsequently found 
useful applications in integrable models without internal symmetry structure. In 
this paper we conjecture the nature of form factor clustering for the general O(n) 
o"-model and make some tests in leading orders of the 1/n expansion and for the 
special cases n = 3, 4. 



1 Introduction 



In this paper we will investigate certain properties of form factors of integrable mod- 
els in two dimensions. This class of models allow a unique field theoretical insight 
because they admit special non-perturbative methods to their solution known as 
the S-matrix bootstrap approach |l|2j . In particular once the spectrum of stable 
(massive) states has been identified a well motivated S-matrix can be postulated. 
Going further one can then attempt to solve equations for the form factors of local 
operators and finally compute correlation functions of these operators by saturat- 
ing with a complete set of intermediate states. There are many examples of such 
applications. Of particular interest are the uncovering of structural relations which 
may have corresponding validity or inspire similar relations in models in higher 
dimensions. 

We will in this paper consider mainly the non-linear O(n) <r-models which have 
the additional interesting feature that they are asymptotically free. In particular 
for the case n = 3 many form factors are explicitly known and one can compute 
vacuum 2-point functions up to rather high energies and compare the results with 
numerical MC and analytic perturbative results [3]. 

One of our main motivations for the present paper arose from our recent work 
on structure functions in these models [3]. A result that particularly intrigues us 
concerns the small Feynman x behavior (Q 2 fixed); we found that this was of the form 
f(x)A(Q 2 ) where the behavior of f(x) at small x reflects the high energy behavior 
of the scattering amplitudes and the function A(Q 2 ) is determined in terms of the 
vacuum 2-point function. We speculated that this structure may be universal in 
asymptotically free field theories, in particular for QCD0. One of the key properties 
in the derivation of the result (for the <r-model) is the property of clustering of form 
factors, and it may be that at least this property has some analogy in QCD. 

Roughly, form factor clustering says that if one considers a multi-particle form 
factor of a local operator and one boosts a subset of particles A uniformly with 
respect to the rest B by a large rapidity A, then to leading order the form factor 
factorizes into a function of A times a product of two functions, one depending 
only on the rapidities of A and the other only on the rapidities of B. The functions 
appearing here are again themselves form factors. 

To our knowledge the first observation of FF clustering was by Smirnov [lj 
for the 0(3) non-linear a-model. Thereafter investigations of this structure were 
mainly pursued for S-matrices without internal symmetry structure (see e.g. [5]- 
|10|). Recently it was used in the construction of local operators in the Sinh-Gordon 
model by Delfino and Niccoli [11]. FF clustering appears to be an extra constraint 

x It holds for example in approximations like naive vector meson dominance. 
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on form factors which can be imposed in addition to the usual axioms. Its main 
application so far has been to identify the operator associated with a particular 
solution of the FF equations i.e. clustering can be useful in model building. 

A systematic study of clustering properties for models with internal symmetry 
has not appeared in the literature so far. It is the purpose of this paper to make steps 
to fill this gap for the case of the O(n) non-linear cr-models. In order to be able to 
make non-trivial tests of the clustering properties of the form factors we first had to 
work out some cr-model form factors explicitly. In particular while considering the 
test in 0(4) we had to work out details of the 3-particle spin form factor and this 
is presented in Appendix D. While the structure (a tensor product of SU(2) form 
factors) was given previously by Smirnov [12], he only considered the case of form 
factors with an even number of particles. We also work out a number of form factors 
in leading orders of the large n expansion both by applying bootstrap techniques for 
this case and by the standard saddle point expansion of the functional integral. A 
by-product of our investigation is the verification of the (expected, but non-trivial) 
equivalence of the two methods. 

The organization of the paper is as follows. In the next section we give a brief 
introduction to the model, in particular the S-matrix is described and form factors 
of some familiar operators are defined. In Sect. 3 we remind the reader of the FF ax- 
ioms, give the 2-particle form factors of operators introduced in Sect. 2 and consider 
the specific case of the 3-particle spin form factor. In Sect. 4, for comparison, we 
first briefly review FF clustering for models without internal symmetry. Thereafter 
we consider the structure of FF clustering in the 0(n) non-linear cr-models. We 
motivate an ansatz for the general form for various cases encountered, and conjec- 
ture a relation of the leading FF clustering behavior to the anomalous dimensions 
of the operators which is reminiscent to the form of operator product expansions 
valid at short distances. In Sect. 5 we consider solutions for the form factors in 
leading orders of the 1/n expansion, and in Sect. 6 we verify that these solutions 
are indeed identical with results derived from the quantum field theoretic formalism. 
Tests of the ansatz in the leading order of the 1/n-expansion are described in Sect. 7 
and tests in the particular cases n = 3, 4 in Sect. 8. Various technical details are 
relegated to appendices. 

2 O(n) non— linear o— model S— matrix and operators 
2.1 The Zamolodchikov S matrix 

Particles in the 0(n) model are characterized by their mass M and the quantum 
numbers (a, 9), where a = 1, 2, . . . , n is an 0(n) vector index and 6 is the rapidity of 
the particle in terms of which the components of its momentum are p° = M cosh 9 
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and p 1 = Msinh#. When two particles scatter there is no particle production and 
the bootstrap S-matrix proposed by Zamolodchikov and Zamolodchikov [13J is of 
the form 



S$(9) = a x {9) 5 cd 5 ab + a 2 {9) 5 c a 5 d b + a 3 (9) 5 d a 5 c b , (2.1) 



where 



and 



°m = 7^ , a 3 (6) = a 2 (9) (2.2) 

ITT — i) 



a 2 (e) = e _ exp{iS{e)} , (2.3) 

where the phase appearing here is given by 

roc 

5(d) = 2 1 — sm(0w)K n {u) . (2.4) 
Jo w 



with kernel 

-TTLU , 

KJu) = — . (2.5) 

We have used the notation x = ^2 * n ^ ne a °ove formulae. 

It is useful to introduce the invariant amplitudes corresponding to s-channel 
"isospin" 1 = 0,2,1: 

S (9) = na 1 (9)+ a 2 (9) + a 3 (9), 

S 2 {9) = a 2 {9) + a 3 (9)=-exp{iS(9)}, (2.6) 
Si{9) = -a 2 (9) + a 3 (6) , 



which obey unitarity Sj(9)Si(—6) = 1. 

The particular cases n = 3, 4 are discussed in more detail in Appendix A. 

2.2 Operators and form factors 

In this subsection we discuss the form factors of the most important operators in the 
model, those of the O(n) spin field, the Noether current and the energy-momentum 
tensor. We also define the form factors of a symmetric, traceless scalar operator. 
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2.2.1 The 0(n) field 

The conventional normalization of the O(n) field is given by its one-particle matrix 
elements: 

(0|$ a (0)|M) = $ ab ■ (2.7) 
The general r-particle matrix elements define its form factors by 

<0|$°(0)|&i, 0i; ...;b r , r r = An /£...&. (01, • • • » *r) , (2.8) 

where 

A 3 = — , A n = l, n>3. (2.9) 

The physical "in" states correspond to the rapidity ordering 9\ > 62 > ■ ■ ■ > 9 r . The 
form factors are originally defined for this ordered set of real rapidities but can be 
extended to the complete complex (multi)-rapidity space by analytic continuation. 
See Sect. El We use the state normalization 

in (a[, 9[; a' r ,e' r \ ai , 9 i; ...;a r , # r ) in = (^) r S a[ai . . . S <ar 6(9[ - 9 1 ) . . . 5{9' r - 9 r ) . 

(2.10) 

2.2.2 The Noether current 

The normalization of the Noether current operators J^ b (x) is fixed by the equal time 
commutation relations 

J a6 (0,x),$ c (0,y) 

where 

t C at = S c a 5 d b -6 d a 5 c b . (2.12) 
The current form factors are given by 

(01^(0)161,01; • • • 'A A)* = ~i^ a q a /ff.^i, . . • A) , (2.13) 

where 

q a = (pi+P2 + ---+PrT, Pi = (p i ,pl) = (Mcoshe i ,Msmh9 i ) (2.14) 

and eoi = — eio = 1- The normalization (|2.1ip implied^ the following result for the 
one-particle expectation value. 

(c,9\J^(0)\d,9) = -2ip^. (2.15) 

2 Recall that for particles with rapidity 6 corresponding to 'bra' vectors in the expectation value 
the form factor functions have to be analytically continued to the complex rapidity value 9 + in. 



= itf d 5(x-y)¥(0,y), (2.11) 
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2.2.3 The energy-momentum tensor 

The energy-momentum tensor is normalized so that its space integral 

/oo 
dxT 00 (0,x) (2.16) 
-oo 

is the Hamiltonian of the system with one-particle eigenvalues given by H\b, 9) = 
M cosh # 1 6, 9). The energy-momentum tensor form factors are 

(011^(0)161, 0i; ...;b r , 9 r )' m = ( M 2 - q^q u ) f bl ... br {0i, ...,0 r ), (2.17) 

where n^ v is the 1+1 dimensional metric characterized by ryoo = — Vn = 1- 
The case n = 3 is discussed in further detail in Appendix B. 

2.2.4 Symmetric, traceless tensor operator 

Finally we define the form factors of a Lorenz scalar and symmetric, traceless iso- 
tensor operator H cd 

(OlE^OJI&i, 0i; . • • ; b r , 0,) in = /^(fli, ...,9 r ). (2.18) 
2.3 Two particle form factors 

Using O(n) symmetry and Poincare invariance, the two-particle form factors can be 
parameterized as follows 

(0\jf(0)\a,a;b,P) = ie^ <? V»i(a - /3)t% , 

(0|S cd (0)|a,a;6,/3) = -i^ 2 ( a - P) t c a f , (2.19) 
(0\T^(0)\a,a;b,P) = ^{q^ - q 2 rj liu )^ {a - (3) 5 ab , 



where 

= % ^ + S d a 6 c b , It = sH-- 5 cd 5 ah . (2.20) 

n 

It can be shown that the normalization of the operators defined above implies the 
following singularity structure for the functions ipi{9) 

—4i 
[9 — m) 1 
2 

faW) » , 9^ in, (2.21) 

9 — ITT 

tp2(9) regular at 9 = m . 
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3 Form factor axioms 



In this section we recall the functional equations pQ satisfied by the scalarized form 
factors, which we generically denote by T a ^...a r (9i, ■ ■ ■ ,9 r ) in this section. It turns 
out to be convenient to introduce the Faddeev-Zamolodchikov operators Z+(9) sat- 
isfying the exchange relation 

z+(9)z+(e>) = slt(o - e')zi{e')z+(0) . (3.1) 

Now we can define the multi-index matrix S& ai ... Or ;& 1 ...& 7 .a(/2|01> ••• ,O r )by the relation 

Z+(f3)Z+(9 1 ) ■ ■ ■ Z+(9 r ) = S bai ... arM ... br a((3\0i,- ■ ■ AK(^) • • • ZtS 6 r)Zt{{3) • 

(3.2) 

The form factor axioms are the following five functional equations [1] 

•7"ai...ar(01> ■ ■ ■ ,0 r ) = jF ai ,,, ar (9\ + A, . . . , 9 r + A) , (3.3) 

T... xy ...{- ■■9,9'---) = S™(9- 9')T... WV ...(- ■■9',0---), (3.4) 

a2-..a r (^1 

+ 2iri, 9 2 , . . . , 9 r ) = F a2 ...a r a 1 (92, ... ,9 r , 6>i) , (3.5) 



lim e T a ba x ...a r (!3 + iir + e,(3,9i,...,0 r ) 

= 2i {$abJ : a 1 ...a r {9l, . . . ,9 r ) — Sb ai ...a r ;6i...fe r a(/3| #1 , • • • , 0r)F bl ...b r (#1 > • • • j ^r)} , 

(3.6) 



•7"ai...a r (#l> ■■■,9 r ) = w p J r a r ...a 1 { — 9 r , ■ ■ ■ , —9\) ■ (3.7) 

In the last equation u; p is the parity of the scalarized form factors. It is equal to 
unity for all operators considered above except for the Noether current, for which it 
is equal to —1. 

Next we define a new type of reduced form factors H by 

T (Q 0\ ^Ol...O r (^l) • • • > ®r) (O o\ 

J-ai...ar{Vl,---,Vr) - n , Q q-s — > 

L> r (pi, . . . , V r ) 

C r (9 u ...,9 r ) = n coshf^-^V (3.9) 

l<i<j<r ^ ' 



3 Here "new" is wrt those usually denned by factoring out the product of 2-particle scalar form 
factors e.g as in Appendix B 
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Three of the form factor equations for F ai ___ ar (6\, . . . ,9 r ) are of the same form as 
(|3.3p . (|3.4p and (|3.7p and the equation corresponding to (|3.5p is only modified by a 
sign factor (— l) r . Finally the residue axiom (13. 6p is rewritten as 

/ ■ \ r r 

F abai ... ar (p + i7r,(3,9 1 ,...,9 r )= TTsinh(/?-^)- 

VV f = \ (3.10) 

{Sba 1 ...a r ;b 1 ...b r a{(3\9l, ■ ■ ■ > #r)-ffei ...6 r (#1 > • • • ; #r) — <5ab-F 0l ...a r (#1 , • • • j #r)} ■ 

3.1 Two particle form factors 

Two of the form factor equations, (|3.3|) and (|3.7|) . are automatically satisfied by 
the ansatz (|2. 19|) . The residue equation (|3.6p does not apply to two-particle form 
factors, but one has to satisfy the normalization conditions (|2.2ip instead. Finally 
flMD and become (7 = 0, 1, 2) 

M8)=Si(6)M-0) (3.H) 

and 

^{0 + 2yri) = (-1) / ^(-9) (3.12) 



respectively where the Si are defined in (|2.6p . 

A useful building block in the construction of form factors is the function 



A(9)=l ^^) COSh[( " + ^ H " 1 - (3-13) 
Jo uj smh 7ro; 

Its main properties are 

A(m + 9) = A(i7r-9), A(in) = Q, A(9) = A(-9) + i5(9) (3.14) 
and its asymptotic behavior for large positive 9 is given by 

A(iw + 9), KeA{9) « —£„(()) - — i^(0) + 0(1) . (3.15) 

2 7T 



If in (|2.4p we substitute K n {uf) by the function k a (io) = — e ™ we get 

e^W = m7T + 9 n . (3.16) 
ia7r — 9 

We denote the related building block by A a (9). 
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Using the building blocks denned above we find 

MO) = tanh^exp{A(#) + A 2x (0)}«ci#-*, 

ih(B) = sinh^ exp{A(6<)} ss c 2 9 x , (3.17) 

MO) = t Z ^iW«c O 0-( 1+ *>. 
Z7T — y 

Here we also indicate the asymptotic behavior of the form factors ipi(9) for large 
positive 9. 

3.2 Three particle form factors 

In this section we write down the form factor equations for the three-particle form 
factors of the O(n) field operators. We start with the definition 

&>(■<*> P>7) = U bV^T^ TT^- ^ 

cosh — ^ cosh ( — tt-M cosh ' 



2 j lAA^ \ 2 ) ^ 2 

The form factor equations (|3.3H3.7j) in this special case become (n > 3) 

^i(«,A7) = i^(a + A,/9 + A,7 + A), (3.19) 
^6c(«,A7) = S^C9-7)i4,(a,7,^), (3-20) 
*lc(a,A7) = F b d ca (P n ,a-2m), (3.21) 

F a t(/3 + iyr, A 7) = \ sinh(/3 - 7 ) {^(/3 - 7) - 5 ab 5 d c } , (3.22) 
= FL ("7, -/?,-«)• (3.23) 

(|3.19p and f)3.21|) are satisfied by the following ansatz 

^o6c( a > A 7) = <5 a Ac^(a - A a - 7) 

+ S bd S ac X(p - 7, - a + 2vri) + 5 cd 5 a6 X( 7 - a + 2irt, 7 - /S + 2m) . (3.24) 

We can also rewrite (|3.22p and (|3.23p in terms of this single function X. We get 

X(u,v) = X(2iri — v, 2ni — u) , 

X(9,in) = ^sinhfl a 3 (9), (3.25) 
X{9,i7r + 9) = -sinhfl [a 2 (9) - 1] . 
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Finally (|3,20p is equivalent to 



X(u, v) = [na\{v — u) + <j2(v — u) + a^(v — u)] X(v, u) 

+ a\{v — u) [X(u — v, 2iri — v) + X(2iri — u, 2iri + v 



(3.26) 



and 



(3.27) 



X(v — u, 2tti — u) = a2(v — u) X(2iri — u, 2ni + v — u) 
+ &3(v — u) X(u — v, 2ni — v) . 

4 Form factor clustering 

4.1 Models without internal symmetry 

We first briefly review clustering properties of form factors for models without inter- 
nal symmetry for which the majority of detailed investigations have been carried out 
so far. In [9] Delfino, Simonetti and Cardy studied models defined as perturbations 
of conformal invariant theories i.e. those formally defined by the action 

A = A CFT +g J d 2 xcf>(x) (4.1) 

where the operator (f>(x) is of dimension 25 < 2. Further attention is restricted to 
those perturbations where an infinite number of integrals of motion survive and the 
resulting massive model is integrable. 

Consider first the case when there is only one species of massive particle of 
mass m and 2-particle S-matrix element S(9). The cluster hypothesis proposes that 
multi-particle form factors of a scaling operator $ of dimension 2<5$ < 2 in such a 
theory factorize according to 

lim Ff +l (6\ + A, . . . , r + A, 6 r+ \, ... , 9 r+l ) = ——Ff(6i, ... , d r )F®(0 r+ i, . . . , 9 r+ i) 

(4.2) 

Actually, (|4.2p is only valid for operators corresponding to primaries in the conformal 
limit. More generally, conformal operators can be classified as 

CmCm$, (4.3) 

where C m is a combination (at level m) of the left Virasoro operators and similarly 
Cm is built from right Virasoro operators. Identifying operators in the massive 
theory with their conformal limit, the generalization of (|4.2[) for descendant operators 
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reads [2] 



lim e - mA F^^{e 1 + A, . . . , e r + a, e r+u e r+l ) 

= F r c ^(e 1 , ...,e r ) Ff-*(e r+1 , e r+l ) 

It follows that knowledge of the multi-particle FF can be used to determine the 
vacuum expectation value (<&) which is in general non-vanishing due to the absence 
of internal symmetries. This is useful since the vev can be obtained by other means 
e.g. by the thermodynamic Bethe ansatz. 

Ref. [9] made the validity of (14. 2|) highly plausible by considering the massless 
limit to the UV critical point in which the mass m —* and rapidities of the particles 
are simultaneously taken to oo such that momenta are fixed. The basic well-accepted 
assumptions are that i) lim^oo S(8) = 1 so that massless left and right movers 
decouple and ii) that the operator space of the conformal point and that of the 
perturbed theory have the same basic structure. In particular to the scaling operator 
<&(x) in the off-critical theory there is an associated conformal operator <1> of the 
same scaling dimension 25$. The property has been noticed in the past to be fulfilled 
by various FF solutions in specific models [5]-[TT], and many examples and tests of 
FF clustering have been successfully performed. 

Tests of the hypothesis for a given S-matrix, involve solving the general func- 
tional equations together with the cluster constraints and seeing whether the number 
of independent solutions equals the number in the corresponding Kac table of the 
associated CFT. Once this has been established one can identify the operators cor- 
responding to the solutions by computing the dimensions <5$ either by studying the 
short distance behavior of the 2-point function computed by saturation by lowest 
states, or using the DSC f')\ sum rule 

6% v - Si R = - JL J d 2 x(Q(xM0))c , (4.5) 

where O(x) is the trace of the energy momentum tensor which is related to the 
perturbing field by O = 4irg(l — 5)4>(x). In a massive theory 6$ R = 0. 

4.2 FF clustering in the non linear sigma model 

To our knowledge FF clustering in the 0(3) non-linear sigma model was first dis- 
cussed by Smirnov p]. A more detailed exposition of this case was presented by 
Balog and Niedermaier [3]. In the following we consider the general 0(n) case which 
exhibits a rather rich structure. 
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4.3 Clustering (leading term) 

Let us divide the particles into two subsets and boost the particles in the first set 
by a large (positive) rapidity A. Form factor clustering means that the scalarized, 
dimensionless form factors have universal large A asymptotics (which usually behave 
as a power in A instead of the constant behavior exhibited by the models considered 
in the previous subsection): 

/a 1 ...o fc 6 1 ...6j(«l + A, . . . ,a k + A, f3i, . . . ,Pi) 

(4.6) 

— hk-i(A) <?oi...a fc ;6i...&j( a l> ■ ■ ■ i a k\P\i ■ ■ ■ j Pi) 
where the clustering function 

M A + *)-M A ) ( 4 - 7 ) 

and also the functional form (and the dependence on A) of the sub-leading terms 
depend on the type of the operator. These sub-leading terms are often suppressed 
by some negative power of A, see Sect. 14.61 

Using the asymptotic properties of the S-matrix, 

S?A0) = 6 a A + ^y i t% + ---, (4-8) 
we can show that the expressions 

5a 1 ...a fc ;6 1 ...fe ; ( Q l» • ■ ■ , a k', Pi, ■ ■ ■ , Pi) (4.9) 



satisfy the form factor axioms (|3.3l )- (|3.6p in the variables corresponding to the first 
set while the dependence on the particles belonging to the second set are playing 
the role of dummy parameters. This is almost trivial for (|3.3p - (|3.5p . and for f|3.6|) 

we get 

faba 1 ...a k b 1 ...b l ('y + ^ + e + A > 7 + A ) a l + A > ■ ■ ■ > a k + A, Pi , . . . , Pi) 
= h k +2;l(A) g a ba 1 ...a k ;b 1 ...b l (l + «7T + e,7>«l> ■ ■ ■ > <*k ] Pi , ■ ■ ■ , Pi) 

2i 

{Sab fa 1 ...a k b 1 ...b l (ai + A, Ojfc + A, Pi,..., Pi) 

~ S ba 1 ...a k b 1 ...b l ;a 1 ..raib 1 ...b l a(. r r + A l«l + A, . . . , a* + A, ft, . . . , fl) 
•4...a fc 6 1 ..i i («l+ A '---«fc + A,/?l,...,A)} 

2i 

— h k -i(A){S ab g ai ... ak - bl ... bl (ai, . . . , a k ; Pi, . . . ,Pi) 

— Sb ai .. (7|«1, . . . ,a k ) ga 1 ...a k ;bi...bii. a U ■ ■ ■ > a kl Pi , ■ ■ ■ , Pi)} ■ 

(4.10) 



e 



11 



From this we see that, because of the uniqueness of the solution of the set of form 
factor axioms, 



h k+2]l {A) = h k;l (A) (4.11) 

and for k > 1 

g a ba 1 ...a k ;b 1 ...b l {'l + VK + e,7,Ol, . . .,Otk,Pl, ■ ■ ■ , Pi) 

2i 



{S a b 9ai.,.a k ;bi...b l { a l, ■ ■ ■ > Pi, ■ ■ ■ j Pi) 

(7|ai, . . • ,ajfe)^a 1 ...a fc ;6 1 ...6j(ai, • • • , afc! • • • , A)} 



(4.12) 

which is nothing but the residue axiom for the first set of particles where the quantum 
numbers of the particles belonging to the second set are dummy parameters. 
Similarly we find that 

h k , l+2 (A) = h k; i{A) (4.13) 

and the axioms (|3.3p — (|3.6|) are satisfied in the variables corresponding to the second 
set, for fixed {ai . . . a k }, {a±, . . . , a k }. 

Using the recursion relations (|4.11j) and (14.131) we see that there are three clus- 
tering families, the cases of odd-odd, even-odd and even-even clustering. 

The function g ai .. Mk bi...bt 1S essentially a product of two scalarized form factors, 
corresponding to the operators B and C. Denoting the original operator by A, the 
clustering relations can be symbolically represented as 

A~h(A)B»C, (4.14) 

or, since in the O(n) model h(A) is always a power, more explicitly as 

A~-Lb»C. (4.15) 

The residue axiom is applicable for clustering for k > 3 only but from [3] we see 
that for O(n) nonsinglets 

MA) = ^ (4-16) 



since 



Fab-M-M (A + in + e, A, ft,... , $) = t% F^ k (fa, . . . , A) . (4.17) 
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where t°^ B are 0(n) generators in the representation under consideration. Similarly 

>ifc;2(A) = 1 (4.18) 



and 



?f 1 ...a fc ;a 6 («i ! ■ ■ ■ , <*kl 7 + «r + e, 7 ) ^ - — /£.. 0fc (ai, . . . ,a fc ) . (4.19) 



For odd, using the uniqueness of the solution, we can solve the problem step 
by step, starting from the k = 1 case: 

ga-M.M^Pi' •••>$) = GW..b ( (/?i, • • • , A) • (4.20) 

Then 

5oi...o fc ;6i...6i(ai, ■ ■ ■ ,0! k ; Pi, . . . , A) = /ai...o fc ( a l) • • • > G abl ,„ bl (Pi, . .. ,A) , (4.21) 

where /^..^ (ai, • • • , «&) is the form factor of the basic field <l> a . 
Similarly for / odd 



9a 1 ...a k ;b 1 -h («1 , • • • , a k \ Pi , ■ ■ ■ , Pi) = #ai...a fc &(ai, • • • , «fc) /b 1 ...6 i (Pi, ■ ■ ■ , Pi) , (4-22) 

where 

0ai...a fc ;&(ai, • • • , »k] P) = H ai _ akb (ai, . . . , a k ) ■ (4.23) 
4.4 Odd odd clustering 

The simplest case is the odd-odd clustering. Starting from 

g a;b (a;P) =T ab (4.24) 

we can build 

fl , oi...a fc ;6i...6i(«l, ■ • • , a k', Pi, ■ ■ ■ , Pi) = T ab fai...a k ( a l> ■■■ , a k) f bl ... bl (Pl, ■ ■ ■ , Pi) ■ 

(4.25) 

In the case of the current operator J cd 

h(A) = ci A - * (4.26) 

and 

Tab - -t c a f. (4.27) 
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Similarly for the symmetric tensor T, cd 

h{A) = c 2 A x (4.28) 

and 

T ab - -K d b . (4.29) 
Finally for the energy-momentum tensor T we have 

h(A) = c A~( 1+x ) (4.30) 

and 

T ab -> - l -S ab . (4.31) 
We can represent the above clustering relations symbolically as 

S ab ~ A x £^$ c «$ d , (4.32) 

T ~ —3— $ a • $ a . 

A 1+ * 



4.5 3 — ► 2 + 1 clustering 

This is the simplest case of even-odd clustering. In the odd-odd case we could afford 
the luxury of using the exact solution of the 2-particle form factors to obtain the 
2 — > 1 + 1 clustering relations. The exact 3-particle form factor is not known for 
general n (except for n = 3, 4 and in the large n limit), but the form factor equations 
can be solved in the clustering limit. We start from the representation 

/: 6c (A,/3, 7 ) = (A - £±2) K - 7 ) + + • • • | • (4-33) 

The functional form of the ansatz (|4.33p is motivated by the following considerations. 
First of all, it is easy to show [5] using the known short distance asymptotics of the 
2-point function and the spectral representation that the spin form factor must not 
grow faster than any power of the momenta, i.e. it is smaller than e eA for any e. 
This observation, the known behavior of the 2-particle form factors and the fact 
that the theory is asymptotically free, taken together make plausible that for large 
momenta also the 3-particle (and higher) form factors vary logarithmically. 
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Now we impose the form factor axioms in the large A limit. This gives for the 
leading coefficient A x abc the equations 

A* 6c (2™ - 9) = A x cb (9) (4.34) 

and 

A x ab M = S v b ?(e)A* uv (-e), (4.35) 

which means that, as expected, A x bc satisfies the 2-particle form factor equations in 
the last two variables. The general solution is 

A* abc (9) = k MO) <$a Sbc + h MO) C + h MO) t%l , (4.36) 

where the kj are constants. In the k = 1 case we also have to satisfy the residue 
axiom 

A* a S* + e) = -^tt c , (4.37) 

which gives ki = —2irx- 

At the next order we find that A x ahc also satisfies (|4.35p but instead of (|4.34|) we 
have in this case 

A^ c (2m -6)- A x acb (9) = -*ttk^ c6 (0) + tori X t% Ke*(0) ■ (4-38) 
Let us define 

KbM = A^ttz -&)- A x acb (9) . (4.39) 

Obviously, 

K x acb {2^-9)+K' x bc {9) = 0. (4-40) 

There is no solution for A x abc unless the right hand side of (|4.38p is also antisymmetric 
under this operation. Thus the consistency of the next-to-leading order gives an 
additional condition on the leading order: 

X tua A vuc + X tuc A vau — k A bac = . (4-41) 

This purely algebraic relation restricts the possible solutions (|4.36p so that only one 
of the coefficients ki can be different from zero and fixes the exponent k as follows: 

1 = case (ko ^ 0) : k = 0, 

1 = 1 case (kij^O): k = 1, (4.42) 

1 = 2 case (k 2 ^ 0) : k = 1 + 2 X ■ 
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Furthermore, the residue conditions give additional restrictions. First of all, since 
ipo (9) has a double pole at 8 = iir, this excludes the ko / (k = 0) solution. In the 
1 = 1 (k = 1) case the residue conditions fix k\ as given above. Finally, since ip2(Q) 
is regular, the coefficient k 2 cannot be determined by this method. 

Putting everything together, we get for the 3^2 + 1 clustering the following 
result: 

2ttv /2tt\ 1+2x 
/* 6c (A, (3, 7 ) - --^ MP ~ 7) C + ixH 2 MP ~ 7) &° + • • • (4.43) 

Here we have reparameterized the constant &?2 (in terms of the new constant H2) for 
later convenience. 

(|4.43p can also be written in terms of the full 2-particle form factors as 

/5, C (A, A 7) = ^ 7) - m f^iP, 7) + • • • (4.44) 

Note that the second piece in the clustering formula (|4.44|) is always subleading 
to the first. Actually, it makes sense only for n > 4, since for n = 3,4 it is not 
dominant over the 1/A 2 correction corresponding to the first term. On the other 
hand, the two terms are close for larger n values and they become degenerate in the 
large n limit. In Sect. 7 we check (14.44j) in the large n expansion. 

Now we calculate these 1/A 2 correction terms in the large rapidity expansion. 
We take k = 1 and get from (|4.33[) 

/* 6c (A, A 7) « ^73+7 { A abdO + ^%bc(Z) + ■■■}, (4-45) 

where £ = (3 — 7. The leading term is: 

Al bc {0 = -2-K X Mi)tt a c . (4.46) 
To calculate the next term we first write 

A x abc (0 = 2zn 2 xMO MZKhc + A 2 (0a?*ac + H0S x Jab} ■ (4.47) 
We also introduce 

A± = A 2 ±A 3 , A = nAi + A+. (4.48) 

of these variables and get 

0+A+(0 = 2, 
C)-A_(0=0, (4.49) 
+ A(0 = 4(l + X ). 



We now rewrite ([4.38P and (|4.35p for in terms 
A-(0 = A_(-{), A_(2jri- 

27T — f 
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From the residue equations we get 

A + (ztt) = 1, A_(«7r) = -2, A(wr) = 2(1 + X ) ■ (4.50) 
This results from expanding the residue equation 

f- bc (a, (3, A) « - S^CS - A)} (4.51) 

C; — ITT 

for large A using the expansion 

S2(A) = « + ^ t%t - ^ Q t™ + ■■■ (4.52) 

It is easy to find the general solution of (|4.49|) and (|4.50p . We get 

A_(f) = ^_(cosh£); £-(-!) = -2. (4.53) 

Since we know that A_ (£) must not grow exponentially for large £ and it should be 
regular we conclude that 

A_(0 = -2. (4.54) 

We also get 

A( ? ) = *^|i -iH±20 (tr + Q. J,_ 1) = (4 . 55) 

ITT — £ 7T 

and requiring regularity here gives 

A(0 = -^ ±2 ^(^ + 0- (4-56) 
Similarly we can find the regular solution of the A+ equations: 

MO = i 2 l lX ~ \ + "(X) cosh | e- A2 *(0 , (4.57) 
tt(1-2x) 2 

where regularity at infinity requires that u>(x) is constant. (|4.57j) is only valid for 
n 7^ 4, because for n = 4 it becomes singular. For n = 4 the solution of the A+ 
equations is 

MO = f (» " {» + J* Q + £) + J* Q - ^) } , (4.58) 
where, again, regularity requires that go cannot depend on the relative rapidity £. 
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uj(x) can be calculated for x = 1 and x = (corresponding to n = 3 and 
n = oo respectively) from the known solutions or for general from the consis- 

tency of the next-to-next-to-leading order equations in the large rapidity expansion. 
Consistency would not fix the value of go for n = 4, and, as we will see later, it is 
actually not a constant, but depends (linearly) on log A. In the light of the l/(n — 4) 
singularity in (|4.57p the presence of the logarithmic term in (|4.58p is not surpris- 
ing. For n = 4 (and only in this case) the 1/A 2 subleading term is accompanied 
by a log(A)/A 2 term. One can see that this is consistent here, because the term 
containing g$ is exactly the same as the one corresponding to the constant H2, the 
coefficient of the k = 1 + 2% term in (j4.43|) . The only way to determine go is to 
solve the full 0(4) form factor equations explicitly. We will consider this problem 
in Section 8 and Appendix D. 



4.6 General clustering 

What we can learn from the three-particle example is that it is useful to include in 
(|4.6|) some of the subleading terms as well. Thus we have to allow for the occurrence 
of several terms (labeled by an index p) of the form 

fai...a k bi...b t (ai+A,...,a fc + A,/3i,...,A) 

^ £ h H ( A ) ^U-M-Jn • • • , «*; A, - , A) , (459) 
p 

where, as before, the coefficient functions g^ have to solve the form factor equations 
for both sets of variable. 

In the case of odd-odd clustering the sum contains only one term so the results 
in Section [4.41 do not change. But the general even-odd clustering formula (k odd, 
I even) consists of two terms: 

J',,''...,,,. /„.../,. ("l + A, . . . , a k + A, . . . , A) 

- fax...a k ■ ■ ■ 1 °>-k) fb™. bl (Pi,..., A) 
•2tt\ 1+2x 



(4.60) 

This can be symbolically represented as 

A A 1+2 x 
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4.7 A conjecture 



We have seen that clustering relations can be represented in the form 

A ~ — B • C (4.62) 

or as a sum of similar terms on the right hand side. We found that the value of the 
exponent n depends on the anomalous dimensions of the operators involved. In all 
cases we studied so far we have the relation 



k = ds + dc — cLa (mod 1) , 



(4.63) 



where 



j 7o 



(4.64) 



Here 70 is the coefficient of the first term (in perturbation theory) of the anoma- 
lous dimension of the operator O and ft) is the coefficient of the first term of the 
perturbative /3-function. 



Explicitly, 



d$ 

dj 
dx 



1 7o = 1 

2 2ft ~ 2 

d T = 0, 

7s 



(i + x) 



(4.65) 



2ft 



= 2 X . 



Here we used the results 



7o 



n — 1 



A) 



n - 2 
4tt ' 



7s 



1 

7T 



(4.66) 



Using this conjecture, it is possible to write down the formula for even-even 
clustering without doing any calculation. For the current operator we get 



fxy 

^ 27TX 

A 



(ai + A, . . . ,a k + A,/3i, . . . , ft 



f xq (n) f yg (8) - f yq 

J ai...a k V~ l J J 6 1 ...&j\r^ / Jai...a k 



09) 



(4.67) 



2tt\ 

a; 



1+4 X 



xq 



(ry) f yq 
a!...a k \ LX JJb 1 ...b l 



09) 



fyq 

J a\...a k ^ 



Here we have used (beyond O(n) symmetry) the fact that isospin form factors 
cannot occur here (they would give a double pole in the residue axiom) and that 
the residue axiom fixes the coefficient of the 1/A term containing the current form 
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factors. The power 1 + 4% is consistent with the conjecture and is uniquely fixed by 
the requirement that the two terms should be degenerate in the large n limit. The 
constant Q is not fixed by these considerations, but by studying the k = I = 2 case 
in the large n limit we can show that 







I +(»( - 

. n 



Analogously for the tensor form factor we get 

(ax + A, . . . ,a k + A, Pi, . 



where 



J ai...a k bi...bi ' 

~ 27T X [ fX 



,pi) 



+ 



A 

27TX 

A 



- /2vr\ 1+2x r~ 



xq 

a\...a k 



-S Xy fZ..a k ^)fZ.M 



n 



n 



-2 + 



(4.68) 



(4.69) 



(4.70) 



5 Bootstrap form factors in leading orders 1/n expansion 

In this section we consider the solution of the form factor equations in leading order 
1/n expansion . We start with the 1/n expansion of the S-matrix and 2-particle 
form factors. 



5.1 S matrix and 1/n expansion 

The 1/n expansion of the coefficients in (|2.ip is of the form 



= o i2 H 1 5- + (J I 

n n z \ 



with 



ai{6) 



2m 



ITT ■ 



a 2 (9) 



2m 
sinh 9 ' 



i = 1,2,3 



a 3 (0) 



2m 



(5.1) 



(5.2) 



bi{6) 



-4-7T / 7T 

i + 



ITT 



sinh# 



-4tT / 7T 

i + 



sinh 9 



(5.3) 



20 



and = T'(z)/T(z)) 



b 2 (9) 



4iri 2n z 
sinh 8 sinh 2 6 



2vr 



Later we will need the combination 

bi(6) + a 2 (6) + a 3 (9) 



+ in ( 2ni 2iri 



ITT 



9 sinh 6 



(5.4) 



(5.5) 



5.2 Large n expansion of the two particle form factors 

The large n expansion of the two-particle form factors is given by the expansion of 
the functions ijjj in (|3.17|) : 



M9) 
M9) 
M9) 



2itanh| f 2vr r , N , „ / 1 
2 J in [ a (6») + 6(0)] + O 



Z7T 



J? 



71" 



tanh 1 <{ 1 + ^ [a(0) + 6(0)] + O ( 



n 



ii 



i <J 1 + — [a(0) - 6(0)] + O ( \ 



n 



n- 



(5.6) 



where 



and 



a{6) 



2?r 27rsinh0 
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1 

47T 



2tt J V 2vr 



2^ 



(5.7) 



(5.8) 



5.3 Large n expansion of the spin 3 particle form factor 



We assume that for large n the function X appearing in (|3.24p has an expansion of 
the form: 



A(n,?j) = 1 ^ h (J ( — 



n- 



rv 



(5.9) 



The form factor equations (|3.25p have to be satisfied order by order in the expansion. 
On the other hand, (|3.26p and (|3.27j) mix the expansion coefficients (beyond leading 
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order). At leading order they lead to 



f(u,v) = [l + oi («-«)] f(v,u), (5.10) 
f(u,v) = f(v,27ri + u) (5.11) 



and at next-to-leading order we have 

g(u, v) = [l + ai(v - u)] g(v, u) 

+ [bi(v - u) + a 2 (v - u) + a 3 (v - u)] f(v,u) (5.12) 
+ a±(v — u) [f(u — v, 2i\i — v ) + f(2wi — u, 2iri + v — u)] 

and 

q(v — u, 2iri — u) = q(2i\i — u, 2ni + v — u) 

(5.13) 

+ a2(v — u) f(2m — u, 2iri + v — u) + 03(7; — u) /(it — u, 2vri — u) . 
5.5.i Leading order solution 

We here summarize the equations the leading order form factor f(u, v) has to satisfy. 

f(u,v) = f(2m-v,2iri-u), (5.14) 
/(0,i7r) = |sinh0a 3 (0) = ^|^, (5.15) 

f(9,iiT + e) = ^ sinh 6 a 2 (0) =tt, (5.16) 

f(u,v) = f(v,2wi + u). (5.18) 
Defining the function 

27T — t> 

with properties 

R(m-e) = 7rsi " hg ; J R(i7r) = vr, R(2m - 0) = R{0) , (5.20) 
we now take the ansatz 

f(u,v) =R(v-u) [s(u,v) + l] (5.21) 
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and verify that (|5.14fj5TT8j) require 

s(u, v) = s(v, u) = s{2ni + u,v) = s(-u, -v), 
s(0,iir) = s(6,in + 9) = 0. 

It is easy to see that regularity and boundedness at infinity allows the trivial solution 
s(u,v) = only leading to the unique leading order solution 

f(u,v) = R(v - u) . (5.23) 
5.3.2 Next-to-leading order solution 

Using the leading order solution (|5.23|) the form factor equations (I5.12p and (15.130 
can be simplified a little. We list here the complete set of next-to-leading order 
(NLO) form factor equations after this simplification. 

g(u,v) = [1 + ai(w - u)] g(v,u) 

+ [bi(v - u) + a 2 (v - u) + a 3 (v - u)] R(u - v) (5.24) 

+ ai(v-u) [R(u) + R(v)] , 

g(u,v) = g(v, 2ni + u) + a2(u)R(v — u) + d3(u)R(v ), (5.25) 

g(u,v) = g(2ni — v, 2tti — u) , (5.26) 

g{6,m) = |smh0& 3 (0), (5.27) 

g{9,iTT + 9) = ^smh6b 2 (e). (5.28) 

We present the NLO solution in several steps in order to make the checking of 
equations (|5.24fj5T28j) easier. We start with 

, , f „, , vk iir IT !T 1 n , . 

g(u, v) = < G(u, v) H — — — — — > R(v — u) 

K ' \ v-u sinh(u-n) R(u) R(v) J v ; (5.29) 

- R(u) - R(v) . 
It is easy to show that (]5.24|) is satisfied if G satisfies 

G{u,v) =G(v,u). (5.30) 

Next we write 

v — u 

G(u, v) = S(u, v) H — r — k(u — v ) — k(v — u) 

sinh(?; — u) 

sinh v { u — in , . . , . , , . ,] 

+ — T7 r { ~. + 2k(u) - k(v + in) - k(v - m) } (5.31) 

smh(!) — u) y%n — v J 

sinh u { v — in 

+ 2k(v) — k(u + in) — k{u — in) 



sinh(u — u) [in — u 
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Here 

*« = 5*(-s)- (5 - 32) 

dOH) and are satisfied if 

S(tt, v) = S(v, u) = S(u + 2ni, v) . (5.33) 
In the next step we write 

in r sinhu (1 + coshu) 



S(u, v) = E(tt, f ) H -7 r \ cosh u — cosh v + 

2sinh(u — u) I smhu 

sinh u(l + cosh v) ~i 



(5.34) 



sinh u J 

In addition to (l5T24"l) and (IST25D . (IST26I) is also satisfied if 

E(w, u) = E(u, u) = E(u + 27ri, u) = E(— u, — v) . (5.35) 

Finally in the last step we represent E as 

„, . cosh it + cosh u / 1\ , , ._ „„. 

E(«, «) = rT ^ R ^ + 3 + *(- 2 ) + «) . (5.36) 

It now follows that all equations (|5,24H5,2"8j) are satisfied if cr(u, v) satisfies the same 
equations as s(u,v) in (|5.22p . 

Although the form factors are bounded and regular functions for all n, the large 
n expansion, as can be seen from (|5.2p . introduces some singularities at rapidity 
differences equal to or in. Nevertheless, one can show that (|5.22p has only trivial 
solution for s(u,v) even if one allows (first order) poles at these special rapidity 
differences. Thus a(u, v) = and the NLO solution is unique. 

5.4 n = oo form factor equations 

In this section we write down the form factor equations in the leading order of the 
large n expansion. In this limit the homogeneous equations take the form 

F ai ...a, r {Qi-> ■ ■ ■ ,0 r ) = F ai „. ar (0i + A, . . . , 9 r + A) , (5.37) 



F... xy -{- ..e,9'...) = - ai (e- e')F... zz ...{- ■ ■ e',e ■ ■ ■ )6 xy + /•:.,,„...{• 



(5.38) 
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.a r (6l + 27TZ, 62, ■ ■ ■ ,9 r ) — ( — l) r 1 ^-OrOl (^2, ■■■ ,9 r ,9\ 



0,10,2- 



(5.39) 



■fox...Or(^l) ■ ■ ■ ,9r) — W p F ar 



(5.40) 



In (15.380 the first term on the right hand side is of 0(1) if the contracted indices 
belong to the same Kronecker delta. Otherwise it is of order 1/n and can be dropped. 

To calculate the residue equation to leading order we first note the recursion 
relation 



Sba 1 ...a r ;b 1 -b r a(l3\0l, ■ ■ ■ , r ) — Sba 1 ...a r - 1 ;b 1 ...b r _ 1 x(l3\0l, . . . , r -l) S x J r {P ~ ® r ) ■ 

(5.41) 

Starting from the r = 1 case 



Sbai;f>ia(/3|#l) = ^aiAifci + ~~ #1 )^afci <5feai 

1 1 / 1 

+ -a 2 (0 - 0i)5 ab 8 aibl + -a 3 (/3 - 0i)5 aai S bbl + O -5 
n n \n z 

it is easy to show by induction that the residue equation takes the form 
nF abai ... ar (/3 + in, (3, 9 U . . . , r ) 



nsinh(/3-^)} { J> 2 (/?-0 fc )}<U Fai...a r {6li ■ ■ ■ ,0r 
3=1 L k=l 

r 

+ ^ Gl (/? - k ) 5 b o k F ai W..a r (0i, ...,9k,...,9 r ) 

k=l 

r 

k=l 



L j • • • ) w k 1 • • ■ > "r J 



Kk 



(5.42) 



(5.43) 



in the leading order of the large n expansion. Again, the last term is of the same 
order as the other terms only if the contracted indices belong to the same Kronecker 
delta and has to be dropped in all other cases. 



5.5 Solution of the leading order equations for the spin field operator 

In this subsection the number of particles, r, is an odd number and we will use the 
notation v = (r — l)/2. For the leading order form factor we take the following 
ansatz 

K...a r (Gl, ...,0 r )=Kj2 5 l W fi 23 • • • K-lr Qr(^(l), ■ • • , a{r) ) , (5.44) 
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where 



Mr 



1 1 



n v 2 v v\ V 2 



(5.45) 



and a runs over the r! permutations of the particles. Finally we used the shorthand 
notation (e is the sign function) 



K 3 = <W (i)M „ R [e (tr(i) - cr(i)) - 6 aU) )] , 
which for physical (real, ordered) rapidities reduces to 



(5.46) 



(5.47) 



We require that the scalar function Q r (9\, . . . , 9 r ) is symmetric under the exchanges 
2 <-> 3, 4 «-> 5,. . . ,r — 1 <-> r and is totally symmetric under permutation of these 
pairs of variables. Further we require that Q r is 2ni periodic in all variables and 
is even under simultaneous sign change of all variables. Then it is almost obvious 
that the ansatz (|5.44p satisfies the homogeneous equations (|O?l45^0lR It is also 
possible to write 



(5.48) 



u(u-l) 



^ai ^a 2 a- A R(&2 — O3) " ' 3a r - 1 a r R\&r-X ~ &r) Qr(&l, ... ,9 r ) + 



where the final dots stand for all similar terms, corresponding to such permutations 
of the variables not leaving the first term invariant. Finally (15. 43ft will also be 
satisfied by (|5.44j) if the set of scalar functions Q r obeys the following three relations 
for r > 3 



Qr+2(0l, P + (3, 2 , . . . , r ) 



n sinh(/3 - 0,) I \J2 



1 



,fc=i 



sinh(/3 — &k 



Q 



r \ai,...,a r 



), 



(5.49) 



Q r+2 (P + iTT,fi,9 1 ,...,e r ) = < 



Y[ sinh(/3 

j=2 



Q r (9i, ... , 9 r ) 



(5.50) 



'Note the relation R(6) = R{-0)[1 + ai(0)]. 
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Qr+2(Pl,P + ™,e 2 ,P,0 3 ,...,0r) 

= sinh(0 2 -9 3 ) I f[ sinh(/3 - 9 S ) \ Q r (9 1} ...,6 r ). 




(5.51) 



Using the results of Subsect. 15.31 we see that 

Q 3 (e 1 ,e 2 ,e 3 ) = i. 



(5.52) 



For r = 5 we have 



Q${01, 02, 03, 04, 05) 



--[l + J2 COsh (0k-0l) 



k<l 



+ cosh(6» 2 + 03-04 



e 5 ) + cosh(0! + e 3 - e 4 - e 5 ) 



(5.53) 



+ COsh(# 2 + 03-01- 05) + COsh(0 2 + 9 3 - 4 " #l) + COsh(0 2 + 01 ~ 9 4 - 65) ■ 

It is easy to check that this satisfies (l5.49fJ53T|) for r = 3. Because of the 2m 
periodicity the scalar function Q r is really a function of the exponential variables 
x k = e^S k = 1, 2, . . . , r. Using the fact that the form factors are regular functions 
that are also regular at infinity and also taking into account the presence of the 
denominator in f|3.8|) we can show that Q r , as function of one of the variables, 
say x 2 , is a finite Laurent polynomial consisting of the terms x^ 1 ^^ 2 , ■ ■ ■ ,x 2 ~ v ■ 
Applying this to the r + 2 case, we see that X2Q r+2 (6i, . . . ,0 r + 2 ) is a polynomial 
of degree r — 1 in x 2 hence it is determined by its values at r different points. If 
Q r is given, these data are provided by (|5.49p and (|5.5ip and we can use them as 
recursion relations to determine Q r + 2 - The solution is given by the explicit formula 



We have applied (|5.54p to determine Qy. We have checked that the function Qi 
constructed this way is also a polynomial in all the other variables and satisfies all 
symmetry requirements together with (|5.50p (which was not used in the construction 
(I5.54p ). It would be interesting to show analogous results for general r. 

5.6 Noether current and symmetric tensor form factors 

The leading order form factors of the Noether current and the symmetric, traceless 
isotensor operator are given by the ansatz 



Kl.a r {0l, ...,0r)=MrY J TaR l4--- K-lr Qr(^(l), • • • , ^(r)) , (5-55) 



Qr+2 (#1 j 02, 03,---, 0r+2) 




i — 1 



(5.54) 



a 
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where r is even, fj,= (r — 2)/2, 

Mr 

and 



1 



2^+V V2 



T ° ~ T a V am a a(2) { d *(l) ~ 0<x(2)) 



where 



i^ a2 sinh# (current), 



3 aia2 



(tensor) 



(5.56) 



(5.57) 



(5.58) 



For even r the scalar function Q r {9\, . . . , 6 r ) is symmetric under the exchanges 
1 <-»• 2, 3 <-»• 4,. . . ,) — 1 r and is totally symmetric under permutation of the last 
/i pairs of variables. Further Q r is 2-7ri antiperiodic in all variables and is invariant 
under simultaneous sign change of all variables. The ansatz (j5.55[) with Q r satisfying 
the above symmetry requirements satisfies the homogeneous equations (15.371 - 157401) . 
We introduce the functions P^ for o = c,t 

Pj c) (#i,... A) (current), 



Q 



r l<71 



, . . . , u r j 



2 cosh 



P®(fii,...,e r ) (tensor). 



(5.59) 



(|5.55|) also satisfies the residue equation (|5.43p if for r > 4 



-isinh(0 3 - 64) < 



H sinh(/?-%) \ P r (°)(0! 



(5.60) 



, . . . , i/ r y , 



and for r > 2 

,(0) 



p% 2 (e 1 ,e 2 ,p + iiT,p,d 3 ,e4,...,e r ) 



1 



-n^-^)|(E smh(/3 _, fc) 



P$°\9 1 ,...,9 rJ . 



(5.61) 



P® 2 (fi + m,O 1 ,p,6> i ,...,0 r ) 

= -i sinh(#i - 2 ) < fl sinh ^ - ^) } P r C)| 



(5.62) 



3=3 
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-2 cosh I — !— - — - ) cosh ^ 



TT , m (5.63) 



j=3 



The functions are given by 



Pi c) {9ue 2 ) = J— P?(9iM = \ (5.64) 

2cosh(^^j 2 



and for r = 4 we have 



T%\KKhM = cosh( gl + g2 ;^- g4 ) , (5.65) 



P?(o u e 2 ,e 3 ,e 4 ) = - cos h ( ) cosh ( ) (5.gg) 

? 2 — #4 



cosh ( — — — ) cosh 



Similarly to the case of field operators discussed in the previous subsection, the 

r-l , . 

expression x 3 2 P^V 2 , treated as a function of the variable X3, is a polynomial of 
degree r — 1 and is determined by its values at r different points. The recursion 
relations (j5.60ff5T63]) can be used to determine this expression at r different points 
and the Laurent polynomials P^ 2 can be calculated from a formula similar to (|5.54p . 

6 1/n expansion of the functional integral 

In this section we check that the bootstrap solutions for the multi-particle form 
factors found in leading order 1/n expansion in Subsections 5.5, 5.6 do indeed cor- 
respond to those obtained by quantum field theoretic calculations [§. 

The 1/n expansion of the functional integral of the O(n) non-linear a-model 
has been described in numerous papers. Starting from bare fields q a one imposes 
the constraint q 2 = n by introducing a Lagrange multiplier field A. Here we just 
recall the resulting Feynman rules for computation of the correlation functions of 
the elementary field: 

q propagator : = 5 ab iD(p,mo) , D(p,itiq) = — ^ , (6.1) 

p z — ?tjq + le 

A propagator : = 2J(p, mo) -1 , (6-2) 

q - A vertex : = -^=5 ab (6.3) 



3 Checks of the S-matrix itself to leading orders in 1/n were performed much earlier [15] . 
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with momentum conservation at each vertex 

- for each external line a factor Z~ 1 / 2 

- for each closed g-loop there is a factor n 

- only g-loops with more than 2 vertices should be drawn 

- integration J over all internal momenta k for which a cutoff A is imposed 
(e.g. Pauli-Villars for the ^-propagator) . 

Renormalization of the bare parameters order by order in 1 jn is given by 

ml = M 2 (l-f;^] ;a s = a s (A/M), (6.4) 



s=l 



oo 



z 



Z = l + Y,^ ;Z S = Z S (A/M). (6.5) 



n 

s=l 



The A inverse propagator function J(q, m) is a special case of the 1-loop inte- 
grals: 

J r (qi,...,q r ,m) = J -^—^ JJ D(k + l j: m) , (6.6) 



where 

Qj = lj — , l-i = lr , (6.7) 

r 

3=1 

which can be computed using the cutting rules. In particular we have 

J(q,m) = J 2 (q,-q,m) = ^—^^ for q 2 = Am 2 cosh 2 (^j . (6.9) 

We will also need the case r = 3: 

hVli, Q2, Qz, m) = — — — * + 2 perms, (6.10) 

D3{qi,q2,q3,m) 

D3(qi,q2,q3,m) = q\q\q\ + m 2 X(q 2 , q%, q 2 ) - ieq x q 2 , (6.11) 

A(xi,x 2 ,x 3 ) = x\ + x\ + x\ - 2x\x 2 - 2x 2 x 3 - 2x 3 xi . (6.12) 



6.1 3 particle spin form factor 

The 3-particle form factor in leading order is, using the Feynman rules above, simply 
obtained from a sum of tree graphs and amputating three of the external lines 
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thereby obtaining: 

/&V 3 (0i>^,03) = ^ 2 * M 2 K^Jfa+Ps)- 1 + 2 perms] + 0(1/ 



n 



(6.13) 



where g = pi +p 2 +P3- Here and in the rest of this section we omit the argument M 
in the functions J r ,D i.e. J(s) = J(s,M). For three incoming (on-shell pf = M 2 ) 
particles one has 

q 2 -M 2 = 8M 2 C 3 (e 1 ,6 2 ,6 3 ), (6.14) 
thus producing Eq. (|5.44p for the case r = 3 with (|5.52p . 

6.2 4 particle current and isotensor form factors 

Consider the current 

J? = q a d»q b - q%q a (6.15) 

whose 2-particle form factor is in leading order just given by the contact diagram 
(q = Pi +P2): 

-^«q a &(On 2) = <bM -P2)» + 0(l/n) , (6.16) 

yielding ffoJfiuO*) = -tf ib2 tanh^l as required. 

The 4-particle current form factor in leading order is a sum of tree diagrams: 

-ie^q v f^ b2b3b4 (0i,O2,O 3 ,e 4 ) = -- $6*6, J (p* + Pi)" 1 

n i<«<i<4 

+0(l/n 2 ), (6.17) 

where k < I and {i,j} U {A;, /} = {1, 2, 3, 4} , q = Ylj=i Pji an d 

V^(q,Pi,Pj) = (2pi - q)^D{q-pi) - (2pj - q) /J ,D(q-p j ) . (6.18) 

For on-shell momenta pi it is clear that q^V fl (q,pi,pj) = as required for current 
conservation, and one can check that 

thus deriving the 4-particle bootstrap solution in subsect. 5.6. 
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Similarly for the 4-particle isotensor form factor: 
cab 



/ftW^fc, ^3, 04) = - £ s tM ~ + ^(9 " 

i<i<i<4 

xJiPk+Phmy 1 + 0{l/n 2 ). 



(6.20) 



For on-shell momenta pi one has: 



D(q- Pi )+D(q- Pj ) 



cosh 



8M^C 4 (6 u e 2 ,9 3 ,9 4 ) 



cosh 



7 fc \ u I ®i ~ 

cosh 



+ cosh ( - — - ] cosh ^ ' 



,(6.21) 



2 J V 2 / V 2 

again consistent with the corresponding result in Subsect. 5.6. 

6.3 5 particle spin field form factor 

The leading contribution to the 5-particle (in state) spin form factor is 0(l/n 2 ). It 
is a little more complicated since there are two types of diagrams contributing: tree 
diagrams with two A-propagators and others involving a closed ^-triangle connected 
to the external lines by three A-propagators. Using the 1/n rules outlined above 
one gets (q = Ej=iPj) : 



fbi 



1 



^ 1 ^ 2 6 3 4 4 6 6 / (5] '(01,02,03,04,05) + 14 perms 



0(l/n 3 ), (6.22) 



where 



/ (5) (0i, 



7 2, (73, C4, P5; 



{£>(pi + p 4 + P5) + D(Pl +P2+ P3)} 
J{P2 +PZ)J{P4+Pb) 



| {D{ P2 + P3 + Pi ) + D(p 2 +P3+ P 5 )} {D(P 2 +PA+ P5) + D(p 3 + p 4 + p 5 )} 



J{q~Pl)J(p2 +P3) 

2 ^(ffl ~ Q,P2 +P3,PA +P5) 

J{q-pi)J{p2 +P3)J{PA +Pb) 
Using (|6.10p this can be rewritten 

r(5)/m Ul 



J(q-Pl)J(P4+P5) 



+ 



u 2 



+ 



u 3 



J(q2)J(Q3) J(qi)J(q 3 ) J(9i)J(?2 



(6.23) 



(6.24) 
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with 



where 



We note 



C/i = Dipt + 92) + D( P1 + g3 ) + D 2 \l 2q3 ] q l , , (6.25) 

-^3191, 92, 93 J 

U 2 = D( P2 + g3 ) + D( P3 + g3 ) + (6 . 26 ) 

Ai(9i,92,93) 

= D(p 4 + g 2 ) + £>(p 5 + g 2 ) + 2(glga)g 3 , (6.27) 

£>3(<7l,<?2,<?3) 



91 = Pi - 9 , 92 = P2 + P3 , 93 = Pi + P5 ■ (6.28) 



D 3 (qi, 92, 93) = 256M 6 C 4 (0 2 , 3 , #4, #5) cosh f^y^J cosh f^y^J , (6.29) 
and then after some algebra it can be shown that 

U 2 = = U 3 , (6.30) 

and© 

JT _ (q 2 -M 2 )Q 5 (9 1 ,e 2 ,9 3 ,e 4 ,9 5 ) 
1 256M4<7 5 (0i,02,03,04,05) ' 1 ' J 

with Q5 defined in Eq. (|5.54|) . thus reproducing the result in Subsect. 5.5. 
6.4 6 particle current and isotensor form factors 

Inspecting the diagrams contributing to the 6-particle current form factor in the 
leading order 1/n expansion we note that we can write these in terms of form factors 
of the spin field (here q = Y^j=\Pj)'- 

v tab 



i£/ivq V fbib 2 b 3 b 4 b B b 6 {01, 02,03, 04, 05, 0&) 

tbic( 2 Pi ~ q)vfb 2 b 3 b 4 b 5 b 6 (02, 03,04, 05, 0a) + 5 similar terms 

tcdiPl +P2 +P3-P4-P5-P6)»fb lb2 b 3 (01,02,03)fb 4 b 5 b (i (04,05,06 



+ 9 similar terms 



+ 0(l/n 3 ). (6.32) 



6 The vanishing of Ui as q goes on-shell is consistent with the absence of particle production. 
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Using the expressions previously obtained we get 



ituvff fbib 2 b 3 b 4 b 5 b 6 (01) 02) #3, 04, #5, ^) = 
(2pi -9)^5(02,03,04,05,06) (2P2 



4^2 { *6l 62 ^>3 b4 ^65 be -ft ( 034 ) # (056 ) 
Q , )m < ?5(01)03)04,05)06) 



C§(62, 03, 04, 05, 06) Cs(5i, #3, #4, #5, 06) 

(gl + P3 + P4 ~ £>2 ~ P5 ~ P&)» _ (P2 +P3+P4-P1-P5- Pehx 

03(61,63, 64)03(62, 65, 6q) 03(62, 63, 64)03(61,65, 6q) 

+ 44 similar terms | + 0(l/n 3 ) . 



(6.33) 



One can check that contracting the rhs with is zero and then obtain the repre- 

(c) 

sentation for the 6-particle current form factor with P 6 given in Appendix C. 
Similarly for the isotensor: 



fb 1 b 2 b 3 b 4 b 5 b 6 (^l^2, 03, 04, 05, 06) 

S6ic/6 2 6 3 6 4 6 S 6 6 (02, 03, 04, 05, 06 ) + 5 similar terms 



+ 



ab re 
s cdJb 1 



b 2 b 3 \ a l> 



, 02, 03)/ 6 l6 5 6 6 (04, 05, 6 ) + 9 similar terms + 0(l/n 3 ) (6.34) 



4~2 { S bib 2 ^>3 b 4 h 5 b 6 ^ ( 034 ) -ft ( 056 ) 

" Q5(02i 03, 04, 05, 06) 05(01) 03, 04, 05, 06) 
-C5(02, 03, 04, 05, 06) C5(01, 03) 04, 05, 06) 
1 1 



+ 



+ 



v C3(0i) 03, 64)03(62, 65, 6q) 03(62,63, 64)03(61,65, 
44 similar terms | + 0(l/n 3 ) . 



(6.35) 



The expression for thus obtained agrees with the corresponding bootstrap so- 
lution as expected. 



7 Large n clustering tests 

7.1 3^2 + 1 clustering in the large n expansion 

We have already computed the 2-particle and 3-particle form factors in the first 
two orders of the large n expansion so we are able to check the clustering formula 
(|4.44p in this limit. 

For the 2-particle form factors we have 



^ I (6) = F I (6) + -G I (6) + ... 1=1,2, 
n 



(7.1) 
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where 



and 



F x {6) = tanh- 



F 2 (9) 



(7.2) 



d (0) = 2vr tanh - [a{9) + b{9)] , 



G 2 {9) = 2m[a(9)-b{9)\ . (7.3) 
In the case of the 3-particle form factors we first write 

Se XU*6 hr .X(A,A + 9) 

(7.4) 



/* 6c (A + /?,/?, 7) 



e y + 



d^oc X(A, 2vri -ff) + 5 x J ab X(9, A + 9)} , 



where 9 = (5 — 7. From here it is easy to calculate that the coefficient of the leading 
(1/n) term on the right hand side of (|T.4[) is 



2vr 
"A 



{F 1 (9)t%-iF 2 (9)s%} + 



Comparing it to (|4.43|) we see that 

H 2 = 1 + X" + 0(x 2 ) 



(7.5) 



(7.6) 



where u is a real number. Using this result and (|4.43|) the predicted form of the 
NLO (1/n 2 ) coefficient on the right hand side of (|7.4p becomes 



2tt 



A 



— U2F 1 + Gx) tit ~ (2F 2 + G 2 - 2F 2 In — + u;F 2 ) + 2,.F 2 ^, 

(Here the argument of all the functions Fj, Gj is #.) 
Now we use the large A asymptotic expansions 



(7.7) 



<?(A,A + 
g(A,2iri- 



ir 

2A l 



l + e e 



7T 



-^e A e* [2 + 2^a{9)] 



+ e y + l 



ta (s) 



2nb{9) 



(7.8) 



5 (0, A + 



-^e A <j [2 + 2vra(0)] 



+ e y + l 



2^6(#)e e 
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and calculate the large A asymptotics of our NLO large n 3-particle form factor 
explicitly. We find that it is indeed exactly of the form (|7.7p , if we choose 



- - 2* [ - 



(7.9) 



7.2 Odd odd clustering 

Specifying the odd-odd clustering formula for the current form factors at the leading 
order in the large n expansion gives 

/„A/2\ kl 

pcd 



K...a kbl ... bl («i + A, • • • , a k + A, ft, . . . , A) 



j.cd 
l ab 



■ cxp < 



I ^ k ^ 



(7.10) 



i=l 



Note the appearance of prefactors which arise because we are considering reduced 
form factors here. 

Analogously for the leading large n form factors of the symmetric tensor we get 

kl 



Fai.a kbl ... bl («i + A, • • • , a k + A, fa, . . . , A) 



3 A/2 



cd 
s ab 



' eXp { oE ' -«Eft' } ^...afcC"!. ■■-.«*) ^...ftjCft.--- »ft)- 



(7.11) 



Using the explicit solution of the current and symmetric tensor as well as the 
spin field form factors we can prove that the k = 1 special case of the clustering 
relations implies 



KKft,ft,...,A) = 

and 

L,(ft,ft>»-.A) 



2 expjifr-ij^i #,,...,#) (7.12) 



1\~ 



exp < 



Qj(ft, ft. •••,#)■ (7.13) 



i=2 



Here the functions Ki and L; are defined by the asymptotic relations 



A(Z - 2) 



2fj(ft,...,A) 



(7.14) 
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and 

P&^A.ft, . . . , A) = ex p{ A( ^ 2 ~ 1) } MA, •••>$)• (7.15) 
We have checked the clustering relations (17. 12() and (|7.13p for / = 3, 5 explicitly. 
7.3 Odd even clustering 

Using the general odd-even clustering formula we can calculate the reduced form 
factor clustering in the large n expansion at leading order: 

F a 1 ...a k b 1 ...b l ( a ^ + A, . . . , Ofc + A, ft, . . . , #) 

27r / e A / 2 \ kl { I k k 1 ] 

~^k\^~) eXP | 2^ Gi " 2^2 Pj [ ( ai '"-' a *) (7.16) 

• { F bi...bi (A, • • • , A) - -^"..6, (A) • • • > A)} • 

(Here A; is odd and I is even.) 

We define the function Mi by the asymptotic formula 

Qi+m,A, fo, . . . , A) = exp | A{1 ~ 2) | M,(/?i,/? 2 ,...,A). (7.17) 

In the special case k = 1 (|7.16j) leads to the following recursion relation: 

1-2 ( l 

e-^M l (l3 1 ,/3 2 ,...,p l ) = (l) 2 exp \ -\ 

W I i=i J (7.18) 

. | 2c osh (^y^) fj (t) 03i,. •• , A) - sinh(/3! - ft) if ...,A) 



We have checked this relation for I = 2, 4, 6. 

Similarly the / = 2 special case corresponds to the recursion relation 

N k {(3i,(5 2 ,ai, ...,a k )= ex P jX^* ~ "T""^ 1 + ^ 2 ' ) | ^fc^ 1 '" • > a *0 > 

(7-19) 



where 

Q fc+2 (A - A, fa - A, ai, . . . , a*) = exp {(A; - 1) A} N k {fa,fa,ax, ...,a k ). (7.20) 
We have checked (|7.19p for = 3, 5. 
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7.4 Even even clustering 

In this case the leading order reduced form factor clustering is of the form 



C.a fc 6 1 ...6,( a l + A '"-. Q! fc + A .^.---.A) 



2vr 
nA 



,A/2 



exp 



1=1 



3=1 



(7.21) 



and 



F oi .a fc 6i...6 ( (ai + A, . . . , a fc + A, /3i, . . . , A) 
/2vr' 



V nA 



a A/2 



exp 



i=i 



+ 2^.. afc («)^.. bi (/3) + 2^(a)^.. 6l G8)} . 
Using the asymptotic relations 

P/^aCai + A, a 2 + A, /3 2 , . . . , A) = exp {(/ - 2)A} X,(ai, a 2 , Pi,..., Pi) 
p/^ 2 (qi + A, A, a 2 + A, ft, . . . , p t ) <= exp | ^ - ^ a| ^(ai,a 2 ,/9i, . . . , Pi) 



(7.22) 



(7.23) 



we have established the recursion relations 



X l (a 1 ,a 2 ,Pi, ...,Pi) 



exp < 



1-2 



(qi + a 2 ) ~y^Pj 



3=3 



2cosh (^) | cosh (^p) P<%, . . . , ft) 



(7.24) 



sinh 



ai — a 2 



sinh 



A -^P/ C )(A,...,A) 
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and 

e 30 J i- a Oyj( Qtlj a 2 , /?!,..., A) 



i-2 
1\~ 



exp < 



i=3 



cosh 



/3i-/3a 



. a x - a 2 \ ( ai - a 2 

sum — 2 cosh 



P«(ft,...,A) ( 7 - 25 ) 



+ ( osh|^) S mh(^)p/ c )(/ 3l! ... ) A) 



We have verified the above relations for I = 2,4. 



8 Clustering in the 0(3), 0(4) models 
8.1 n = 3 

We recall the discussion in [3]. One has for the reduced form factors (see Ap- 
pendix B): 



9b l ...b m a 1 ...a k (Ai> • • • > A», «i + A, . . . , a k + A) 



A 



km—l 



eabc9° ai ...aM> ■ ■ ■ ><*k)9b 1 ...b m (Pi> ■ ■ ■ > A») + 0(A km ~ 2 ) , (8.1) 



and similarly 



/ifei...6 m ai...a fc (/3l, • • • , An,«l + A, . . . , £*fc + A) 

f/3 1 ,...,/3 m ) + 0(A fem - 3 ). 



a km— 2 a 
9a\. 



.a 



ai, . 



>®k)g bl ...br, 



1.2) 



Note that in (|8.1|) members of the isospin 1 family are mapped into themselves, 
while in (|8.2|) members of the isospin 1 family are linked to members of the isospin 
family. Observe also that there is no distinction between the factorization properties 
between even and odd members of the same family. 

For the special case of k = 1 , m = r — 1 the clustering relations read 



<&...«,. (ft,..., A) 

h ai ...a r (Pl, ••• , AO 



r 2 wi..^(ft, • • • , a-i) + our 3 ) , (8.3) 

ir—3„a r (a Q \ \ r\l qt—A\ 



0r 9t..ar-APl,- + 0(fi-l 



(8.4) 



which is in accordance with the propery that the reduced form factors are polynomi- 
als of partial degree (r — 2) and (r — 3) in the isospin 1 and case, respectively. Since 
the product *$> r also factorizes under clustering, the full (scalarized) form factors also 
satisfy clustering relations, which are similar to (18.11) and (18. 21) . For the I = 1 family 
they can be found in Smirnov's book pp. 



39 



The clustering relations closely resemble some classical equations satisfied by the 
operators. For example dividing an even number of particles into two odd clusters, 
(|8.ip can be interpreted as the quantum counterpart of the current in terms of the 
spin operators. The division of an even number of particles into two even clusters on 
the other hand, resembles the classical equation d^J^ — dyJ 1 ^ oc e a b c J^J£, . Finally the 
clustering of an odd number of particles corresponds to d^s a oc e a b c s b J^. Similarly 
(|8.2p corresponds to the defining equation for the energy momentum tensor in terms 
of the spin fields or equivalently to its Sugawara form T^ v oc J® J" — ^rj^JpJp. 

In ref. [3] the clustering properties above were used in particular to deduce the 
clustering properties of absolute squares of the form factors, summed over internal 
symmetry indices which enter in the expressions for spectral densities. 

8.2 0(4) form factors: a three particle example 

Just as the 0(4) S-matrix (jA.19|) is (minus) the tensor product of two chiral Gross- 
Neveu S-matrices, the 0(4) form factors can be written as tensor products of two 
chiral Gross-Neveu form factors. More precisely, the 0(4) form factors can be 
written as linear combinations of several such tensor products. This solution of the 
0(4) form factor equations, for the case of even particle numbers, was given by 
F. A. Smirnov [12]. The odd particle form factors must have a similar structure. 
The solution for the three-particle form factors of the 0(4) field operator was found 
by M. Karowski [16J: 

fp i ABc(0 1 ,6 2 ,0 a ) = D(p li 2 ,es)T, ^lih^hMF^^jOuM). (8.5) 

UJ 

Here ^y^, c (0i, 02, #3) for w = ± are the spin s = ±| SU(2)-symmetric chiral Gross- 
Neveu model form factors discussed in Appendix D. They satisfy the following ho- 
mogeneous bootstrap equations: 

F^ t2l3 (e 1 + 27ri,e 2 ,e 3 ) = -iw^fefljA) (8.8) 

and the residue equations 

^ 03) - a Sp_ m + ^ SfkV - ^(03)} • 

(8.9) 

Here the chiral Gross-Neveu S-matrix Sfj(9), the anti-symmetric charge conjugation 
matrix and the one-particle form factors F^ (9) are all defined in Appendix D. 
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It is easy to show that the homogeneous form factor equations are satisfied by 
(|8.5p if the scalar prefactor D is shift-invariant, anti-symmetric under the exchange 
of any pair of rapidities and is 27ri-periodic in all rapidity variables. It also has to 
have a first order zero at points where two rapidities differ by in in order to satisfy 
the residue equation as well. The solution is 

^i^ 2 ,03)=9n coth (^) > ( 8 - iq ) 

i<j ^ ' 

which also has the right normalization. With this choice we have 

fp-ABc(a,P,0) _ - _ . ^2 \c aibl c PlCl + iuc aikl Cp lh S^ ll (l3-9)\ 

"| C a2fe2 C P2C2 + C a2fc2 C P2h ^b2C2 (@ ~ ^) } 

= a — ~\ Cab Cpc - C ak Cpl S bc ((3 - 9) \ . 

a — p — in y y 



.11) 



Here the 0(4) charge conjugation matrix Cab is defined in (jD.lOj) . 
We parameterize the 0(4) form factors as follows: 

fp;ABc(@i, #2, #3) = CpaCbc 9i(Q\, 62,63) 

+ C PB Cac g2(Oi,02,9 3 ) + Cpc Cab 93 (61 , 02 ,03 ) • 

Using the tensor product solution we can write this as 

g 2 +g 3 = 2DF { _ +) F { _- ] , 

92-93 = D + F W F[- ] ) , (8.12) 



A 9l +g 2 +g z = 2DF^ ] F { + >. 

For (0i, 02,0^) = (A, a, (3) and using the large A expansion of the components F±^ 
described in Appendix D we find for the 0(4) form factors 



* ( - + -L\ + $ (I _ JL 
2m) \2 2m 



92+93- ~^( i7T ~ 0^1(0 1 9o + ~ 

2^1(0 f 1 . in . 
g2 ~ g3 ^~ A -^ I A 

4ffi +92 + 93 ~ ^2^1 (O(C + «0 + ••• > 

where £ = a — /?. We see that this is exactly the same expansion as the one we found 
by directly solving the 0(4) form factor equations asymptotically in Subsect. 4.5. 
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Only the value of go cannot be determined from the asymptotic solution. This we 
found in Appendix D by expanding the complete solution for large A: 



9 Concluding remarks 

In the course of this work various intriguing relations concerning form factor clus- 
tering in the O(n) sigma-models were discussed and new structures revealed. The 
relationship of the pattern of clustering to the classical field equations in the case 
of 0(3) has been previously known |l|3j . Some of these patterns in particular those 
involving the Sugawara structure of the energy momentum tensor and those involv- 
ing the (non-Abelian) curl-freeness of the Noether current (which is so important 
for integr ability) probably extend to general n. We have further formulated a con- 
jecture in Subsect. 4.7 concerning the (on-shell) nature of clustering to the operator 
product expansion. 

We have tested our ansatz in various examples. Firstly we checked that the 
solutions obtained by solving the form factor equations in leading order 1/n coincided 
with those obtained by the field theoretical approach to the model. Although this is 
as generally expected, it constitutes yet another test of the proposed equivalence of 
the S-matrix bootstrap construction and functional integral definition of the models. 
We also found that the large n limit and limit of large rapidity commute. Although 
this is observed in previous studies it is not an obvious fact (recall that there are 
many examples where the large n limit and limit of small rapidity do not commute). 

The case n = 4 is a special case. Here we studied in detail (to our knowledge 
for the first time) the 3-particle spin form factor. The tensor product structure is 
probably particular to this case but its general form involving the hypergeometric 
functions may give a hint to the outstanding unsolved problem of the construction 
of form factors for general n > 4. Moreover in this case we found that subleading 
terms in the form factor clustering involved also logarithms of the (large) rapidity 



As mentioned in the introduction it is not completely implausible that some of 
the structural properties found concerning form factor clustering in integrable 2-d 
asymptotically free models have their analog in 4-dimensional models in particular 
for processes when the kinematics effectively reduces the dimension to 1+1. 
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Appendix A. Explicit S— matrices 
A.l The n = 3 case 

For n = 3 the kernel K$(uj) = e _7ra; and the integral in f|2.3[) is easily done. The 
result is well known [13]: 

a x {0) 

2 (0) = ll} ,'";i, /fJ ' J ' (A.l.) 
0-3(0) 

A. 2 The n = 4 case 
In this case (|2.2|) and (|2,3|) simplify to 



2m9 


(6 + iir)(6 - 


2m) 


9(6 - iir 


) 


(6 + iir)(0- 


2m) 


2iri(i7r — 


0) 


(0 + m){6- 


2m) 



2 (0) = ^—S {2 \e), SW(8) = -A 2 (6), (A.2) 

— ITT 

03(0) = -^S (2) W, 



^7^ 



where 

The 0(4) S~-matrix is here given in the real basis 

\a,0), o = l,2,3,4. (A.4) 

It is useful to transform the particles into a complex SU(2) x SU(2) basis 

\A,0), A = ++,— ,+-- + . (A.5) 

The transformation and its inverse are given by 

\a,9) = n a A \A,6), \A,6) = lC A a \a,8), (A.6) 
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with /c A «n a B = «5f , n a A )c A b = s b . 

The transformation rule of the S-matrix is 

s%&0) = ic A 'ic B b n e c n i D s3{0). 

Using the 0(4) S-matrix explicitly we get 

S C ab(0) = <Ti (0) P CD Rab + a 2 (9) 5 C A 6% + a 3 (9) 5 C B 5% , 

where 



P CD = n c n D 



X " U X 1 



Rab = K-a x &i 



(A.7) 
(A.8) 
(A.9) 



Next we define the SU(2) x SU(2) basis explicitly. The 0(4) generators in the 
vector representation are 



We now define 
1 
2 

Further 



fr ab ^j = i [d ax 5 by - 5 ay 5 bx ^j , 
and A k = t m , 



a, b = 1,2,3,4. 



k, I, m = 1, 2, 3 . 



These are the SU(2) x SU(2) generators since 



wLwl 



o, 



We now define SU(2) x SU(2) particle states such that 

| + +, 9) has eigenvalue \- w.r.t W\ and i w.r.t w£ , 

\-\ — ,6*) has eigenvalue ^ w.r.t W\ and — ^ w.r.t 



(A.10) 

(A.11) 
(A.12) 

(A.13) 

(A.14) 
(A.15) 



and so on. Here the phases are also important. We make the following choice 

| + +,0) = -L{|l,0)-i|2,0>}, 

\--,9) = ±{\l,9) + i\2,e)} , 

\ + -9) = -j={i\3,9) + \4,0)} , 

\-+,e) = -L{i\3,9)-\4,9)} . 
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Now we calculate 

P cd =oj c 5 cD (no sum), R AB = ua5 A b (no sum), (A. 16) 

where 

lo ++ = uj = 1 and lu^ = u (. = — 1 (A. 17) 

and B is the charge conjugate of B. 

A. 3 Tensor product S matrix 

The S-matrix of the SU(2) chiral Gross-Neveu model is [1] 




(A.18) 



where a, 0, 7, <5 = +,— . Taking the tensor product of two such S-matrices we get 
~Sl\i(0) Sll%(9) = a 2 {6) gn 8% 5% + a 3 (6) 8}\ dg 5 s a \ 8% 

+ °m (<$ * - ( 5 l 4 - *2 • (A ' 1! 

which is the same as (|A.8j> if we make the identification A = (01,02) etc. 
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Appendix B. Some form factors for the case n = 3 

For n = 3 we can define J® = \e abc J bc and rewrite (|2,13p as 

<O|J2(O)|6i,ei;...;^,0 r ) ta = -te pa 9 a, /£.. Ar (6> 1 ,...,(9 r ). (B.l) 

Using this unified notation^ for the form factors of the 0(3) field and the Noether 
current we can introduce reduced form factors g% br by 

&M 1> • • • A) = *r(01, • • • ,0r)g a bl ... br (Ou ■ ■ ■ , O r ) , (B.2) 



where 



with 



* r (0 1 ,...,O r ) = ^ir&- 1 ) J] (B.3) 

l<i<j'<r 



Similarly for the energy-momentum tensor we define 

f bl ...b r (01, ...,0r) = Vrtfl, • • • A) g bl ...bM, ■ • • A) • (B.5) 

The advantage of using these reduced form factors is that they, with only one ex- 
ception, are polynomial expressions in the particle rapidities [3]. 

The first reduced form factors for the 0(3) field are g% {Oi) = 5 abl and 

+ 5 ab2 5 blb3 (6i - 3 - 2iri) + 5 abl 5 b2ba (d 3 - 2 ) . ' ' 

For the current we have 

<6 2 (#lA)=6^ (B.7) 

and finally for the energy-momentum tensojfl 

9b lb2 (e l ,e 2 ) = - — - — -. (b.8) 

"l — "2 — 27T 

Note that for n = 3 the functions ipo,tpi defined in (|2.19p are given by 

^l(fl) = - T m ^0(0) = ^^. (B.9) 

Many further explicit examples can be found in ref. [3]. 



7 Note that no confusion can arise here since the O(n) spin field has non-vanishing form factors 
for odd number of particles only whereas the current form factors are non-vanishing for an even 
number of particles only. 

8 This is the only exceptional, non-polynomial reduced form factor. 
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Appendix C. Current 6-particle function 

^ 6 {c ViA,#3AAA) = ^[{ 





3coshi(0i + 2 + 3 - 

2 


-04- 


05- 


6 )+3coshi(0i+ 2 + 05 - 3 
Z 


-04 


" 06) 


+ 


3 cosh - (6»i + (9 3 + 4 - 
2 


-02- 


05- 


6 ) + 4coshi(0i+ 3 + 05 - 2 
2 


-0 4 


" 06) 


+ 


2 cosh 1 (0i + 30 3 - 2 
2 


-04" 


-05- 


-0 6 ) + 2cosh^(0 3 + 30i-02- 


04" 


05" 


06) 


+ 


2 cosh - (0i + 30 5 - 2 
2 


-03- 


-04- 


-0 6 ) + 2coshi(0 5 + 30i-02- 
2 


03" 


04- 


06) 


+ 


2 cosh i (0 3 + 30 5 - 0i 
2 


-02" 


-04- 


-0 6 ) + 2 cosh I (0 5 + 30 3 - 01 - 
2 


02" 


04- 


06) 


+ 


2 cosh - (0 3 + 30 4 - 0i 
2 


-02" 


-05- 


-0 6 ) + 2 cosh - (0 5 + 30 6 - 0i - 
2 


02" 


3 - 


04) 


+ 


2 cosh -(0i + 30 2 - 3 
2 


-04- 


-05- 


-06) 








+ 


2 cosh - (0i + 2 + 30 3 
2 


-04- 


-05- 


- 30 6 ) + 2 cosh - (0i + 2 + 30 6 - 
2 


-04" 


-05" 


- 30 3 ) 


+ 


4 cosh - (0i + 30 2 + 3 

2 


-04- 


-05- 


- 30 6 ) + 2 cosh - (0i + 302 + 3 - 
2 


-05" 


-06" 


- 30 4 ) 


+ 


4 cosh -(0i + 30 2 + 5 
2 


-0 3 - 


-06- 


- 30 4 ) + 2 cosh - (0i + 30 2 + 5 - 
2 


-03" 


-04" 


- 30 6 ) 


+ 


4 cosh - (0i + 3 + 30 4 
2 


-02- 


-05- 


- 30 6 ) + 2 cosh I (0i + 3 + 30 4 - 
2 


-302 


"05 


-06) 


+ 


2 cosh - (0i + 5 + 30 6 


-30 2 


-0 3 


-0a) 








+ 


cosh - (30i + 30 2 - 3 


-04- 


-05- 


- 30 6 ) + cosh - (30i + 302 - 3 - 


05" 


06" 


30 4 ) 


+ 


cosh i (30 3 + 30 4 - 0i 


-02" 


-05- 


- 30 6 ) + cosh i (30 3 + 30 4 - 01 - 


05" 


06" 


30 2 ) 


+ 


cosh - (30 5 + 30 6 - 0i 


-02" 


-0 3 - 


- 30 4 ) + cosh ]- (30 5 + 30 6 - 0i - 


3 - 


04" 


30 2 ) 


+ 


cosh - (0i + 30 3 + 30 4 


-02" 


-30 5 


- 30 6 ) + cosh 1 (0 3 + 30i + 30 2 


-04- 


-30 5 


- 30 6 ) 


+ 


cosh - (0 5 + 30i + 30 2 


-06" 


-30 3 


- 30 4 )} 









+ 



{03 ^ 04)} + {05 ~ 06)} + {03 " 04 , 05 ~ 6 )} 



+ 



(C.l) 
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Appendix D. XY model and chiral Gross— Neveu model form factors 

The bootstrap solution of the XY model [T7] is based on the extremal Sine-Gordon 
S-matrix, which is given by 

SZZ(9) = S++(9)=A(9), (D.l) 
S+Z(9) = SZ+(6) = -— e A(8), (D.2) 

Stl(9) = S+Z(9) = -^-A(9), (D.3) 

ITT — 9 

where A(9) is given in (|A.3j) and k = 1 for the XY model S-matrix. The choice 
k = —1 gives the SU(2) symmetric chiral Gross-Neveu S-matrix (|A.18|) . We will 
keep the notation S^L(9) for the XY model S-matrix and will denote the chiral 

Gross-Neveu S-matrix by S^(9). 

For the XY model the crossing relation is 

Sl & ^-9) = c^c &v SZ{9), (D.4) 
where the charge conjugation matrix c a p has non-vanishing components 

c+_ = c_+ = 1 , (D.5) 
whereas for the chiral Gross-Neveu case crossing is given by 

S^(i7:-9) = c^csuSZ(0), (D.6) 

with 

c-| = — c |_ = i . (D.7) 

As discussed in Appendix A, the 0(4) S-matrix is (minus) the tensor product of two 
chiral Gross-Neveu S-matrices: 

S%E(0) = -Sl[%(9) Sl\%(9) , (D.8) 

where A = (a.\,a.2) etc. The crossing relation is 

S%B(i* ~9) = C BM C dn S c a ^{9) , (D.9) 

with 

Cab = c ai p 1 c a2 f3 2 ■ (D.10) 
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D.l SU_i(2) symmetry 

As is well known, the Sine-Gordon model has quantum group symmetry SU g (2), 
which becomes SU_i(2) in the extremal case. In this case the algebra of generators 
t±, j is identical to the ordinary SU(2) algebra: 

[T+,T-] = 2j, U,T±]=±T ± , (D.ll) 

it is only the co-product A that is different from the classical case: 

= j <8> 1 + 1 <8> j, A(r±) = r± <g> (-l) 2j + 1 <8> r± . (D.12) 

This means that if we build tensor product representations from the basic doublet 
representation, the representation matrices are given by 

A 2 (t±) = -t± (g> 1 + 1 (g> r± , (D.13) 



A 3 (t±) = t±®101-10t ± ®1 + 181®t± (D.14) 

etc. These are representation matrices of the classical SU(2) algebra and they are 
simply related to the usual ones. For the two-particle case the relation is 

I + +)cl = I + +), 

l + " >cl = l + " ) ' (D.15) 

I -+)d = -!-+), 

where the \aj3) c \ states transform according to the usual two-particle representation. 
Similarly for higher states we have 

r 

\a r ... a 2 a 1 ) c i = {ai) l+1 \a r . . . a 2 a{) . (D.16) 
i=i 

We denote the SU(2) generators acting in the Hilbert space by J, T±. We will 
be looking for local fields 4>±(z) transforming as elements of an SU(2) doublet: 

2[j,<f> ± (z)\ = ±<f> ± (z), [f + ,Mz)\ = lf-,cp-(z)] =0, 

[f + ,4>-{z)\ =4> + {z), [f-,M*)] = 4>-(z). 

Note that ((f>+)^ / op- for any constant c (the above equations cannot have such 
solutions) i.e. the doublet fields must be genuinely complex. 
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SU(2) symmetry restricts the form factors of doublet operators. For the 1— 
particle form factors we have 

F p , a (9) = ((%(0)|a, 9) = iG(9)c pa . (D.18) 

For the 3-particle ones 

F p - abc (^ d 2,0 3 ) = (0\cj> p (0)\a,e 1 ;b,d2;c,e 3 ) (D.19) 
we introduce the notation 

— (#1> S2, #3) = Fi(6i,02, 9 3 ) , 
F-. + - + (9 1 ,02,0 3 ) = F 2 (0 1 ,92,e 3 ) 1 (D.20) 

-F-;-++(01> #2> ^3) = F3{6l,92, 9 3 ) . 

Note that all other components either vanish by charge conservation or are related 
to these ones by charge conjugation: 

Fp;abc(Ql>®2, O3) = — F p . abc (9i, 02, O3) . (D.21) 

The restriction coming from SU(2) symmetry is 

F 1 -F 2 + F 3 = 0. (D.22) 

We also introduce the form factors corresponding to the manifestly SU(2) in- 
variant basis (|D.16j) : 

F p -M0i,9 2 ,e 3 ) = (O\cP p (O)\a,9 1 ;b,02;c,0 3 ) cl . (D.23) 

For these form factors we have 

Fi = F u F 2 = -F 2 , Fs = F 3 , (D.24) 

F mld (9i,9 2 ,e 3 ) = F p . abc (0i,9 2 ,9 3 ) (D.25) 

and 

Fx + F 2 + F 3 = . (D.26) 

We note that the basic spin fields [17] of the XY model, 

S ± (z) = S 1 (z)±S 2 (z) (D.27) 

obviously satisfy (5 ,+ )^ = S~ and hence cannot be elements of a doublet. 

In the following we will consider the form factors of not only the doublet opera- 
tors but also more general, charge —1, spin s fields. We will use the notation (|D.20|) 
also for these more general form factors. Of course, (|D.22p only holds for the SU(2) 
doublet case. 
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D.2 3 particle form factor equations for general spin 

We recall the bootstrap equations satisfied by the form factors of a charge —1, spin 
s operator (which may or may not be the lower component of an SU(2) doublet). 
The homogenous equations are: 

Fi 1 i 2 i 3 (9i + X,92 + X,9 3 + X) = e sX Fi 1 i 2 i 3 (6i, 62, 9 3 ), (D.28) 
Fiii2i 3 (8ii 62, 63) = S™i A (P2 — 9 3 ) Fi lVU (0i, 9 3 , 62), (D.29) 
F ili2i3 (9 1 + 2Tri,0 2 ,9 3 ) = % i^^sA) ■ (D.30) 

This is supplemented by the residue equation 

F hi2i3 (a,(3,0 3 ) « - _ j- - [c ill2 F l3 (9 3 ) - Vn c ilk Sg^tf - 9 3 ) ^3)} • (D.31) 

Here and in the cyclic equation (|D,30P rj^ is a phase factor that expresses the rela- 
tive non-locality between the field, whose form factors we are constructing and the 
basic spin fields that create the asymptotic particles using the LSZ asymptotic 
formula. Consistency between the cyclic equation (|D.30P and the shift equation 
(|D.28|) requires that 

V+ = r, = e 2wis , r?_ = e~ 27ris (D.32) 

and this is sufficient to determine the one-point function 

F j (9)=g5+e se (D.33) 

up to the normalization constant g. 

If we write the shift and cyclic equations in terms of the independent compo- 
nents Fx, F2, F 3 we get 

F k (9 1 + X,92 + X,9 3 + X)=e sX F k (9 1 ,9 2 ,9 3 ), ft = 1,2, 3 (D.34) 

and 

F k (9 1 + 27ri,02,9 3 )=r ] F k+1 (02,9 3 ,9 1 ), k = 1,2. (D.35) 

The k = 3 equation ^3(6*1 + 2iri, 92, 9 3 ) = tj~ 1 F\{92, 9 3 , 9\) is already a consequence 
of the above two. 

We have seen that SU(2) symmetry requires 

( = F 1 -F 2 + F 3 = 0. (D.36) 

For later purposes we introduce 

F ± = F 1 ±F 2 , (D.37) 
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in terms of which (|D.36P can also be written as F 3 = —F^. From (|D.35|) it follows 
that 

C(o 1 + 27ri,9 2 ,e 3 ) = -TiC(O2,o 3 ,0 1 )+(v + ^\ F 1 (9 2 ,e 3 ,e l ) (d.38) 

thus an SU(2) doublet field must have r/ = ±i i.e. spin s = ±1/4 (mod 1). We will 
consider two such doublet solutions <f>^\z) with s = u>/4, 77 = ui (u = ±). It is 
natural to write the form factors in this case using the manifestly symmetric basis 
vectors (|D.16p , We choose the normalization (g^ = 2) such that the 1-particle 
form factors are given by 

= (0\^H0)\a,d) = 2ic pa e^/ 4 . (D.39) 

Written in terms of the form factors of these operators, the three-particle equa- 
tions (|D.28llD.31j) become the form factor equations (18, 6118. 9p . discussed in the main 
text. 



D.3 Reduced form factors 

To simplify the solution of the 3-particle form factor equations we introduce a set 
of "reduced" form factors f m (m = 1, 2, 3) by writing [17] 



F m (e 1 ,e 2 ,e 3 ) = -27r 2 MY(e 1 ,6 2 ,e 3 )e s ^ +M ^ f m {9i,e 2 ,e 3 ) , (d.4o) 

where the prefactor Y is 

Y(e 1 ,9 2 ,9 3 ) = l[y(9 i -9 j ) (D.41) 

i<j 



with 



where 



,(9) = sinh (0 e E ^ , (D.42) 



= do; [ooA^Tr + ig)-!] 1 ( } 

' uj sinh7rw (l + e™) 



Finally the normalization constant is 



N = —^e- EW e- i7TS . (D.44) 

7T ' 
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Note that the function E(9) is related to ipi(0) (defined in (|3.17p ) for the case n = 4 
by 



V>l(#)|n=4 = 




(D.45) 



and for large it behaves as 

AE(6) = -0 + ln(20) + 2E(0) + m (l - 0~ l ) + O(0~ 2 ) . (D.46) 
For later use we introduce the function 

*^= r G + ^) r (-^)- iDA7) 

With the help of this function we can write the S-matrix element A(0) as 

m = ^iy (D.48) 

For completeness, we give here the form factor equations, rewritten in terms of 
the reduced form factors f m : 

f m (0i + A, 2 + A, 6*3 + A) = e" 2sA f m (e u 2 , 3 ) , (D.49) 



f 3 (a,e,0') = f 3 (a,e',0), 
f-(a,e,0 , ) = U(a,e',0), 
i-K + 9 



f± = h± h 



(D.50) 



f + (a,6,0') 



17T 



7 U(a,0',9) 



fm(01+27ri,0 2 ,03)=fm+l(02,03,0l), ("1 = 1,2), (D.51) 



A (a, (3, 0) » ^- — <^ — - $ /3 - 0) 

a — p — 27T [ 27T — p + 
a — p — Z7T [ Z7T — p + ' 

/s(a, /3, 0) ~ % . <&(0 - /3) - ~ <&(/? - 6 

a — p — Z7T [ 77 

Later we will explicitly solve the form factor equations for the reduced form 
factors and calculate their large rapidity limit. But even before having the complete 
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solution, a lot of information about their large rapidity behavior can already be 
obtained by expanding the equations themselves. We take the Ansatz 



/m(A,a,/3) 



e^V" 
(A - o)p 



W m (a-f3) , 

U m {a - f3) H h 

A — a 



(m=l,2,3) 

(D.53) 



for large A, where k = \ — 2s. We get restrictions on the leading power p and the 
expansion coefficients U m , W m by substituting this Ansatz into the reduced form 
factor equations. In particular, from the residue equations, using the asymptotic 
expansion of $, we get for large A 



f m (a,(3,A) 



a — (3 — iir 



V 



3i7T7] 



(A - /?) 3 /2 4A 5 /2 
1 + ir] 1 



+ 



vr (A-/?)V2 



3n — i 1 
j • 1_ . 

4 A 3 / 2 



- I 1 - - 

7T 



r?y (A-/3)V2 4 V?? 



1 



+ 



where 



c = 4ir v 27re 4 g . 



A 3 / 2 

(D.54) 



(D.55) 



From here we see that for most spin values for which r] ^ i (—\ < s < \), the leading 
power is p = 1/2 but for s = \ we have rf = % and the leading power is p = 3/2. It is 
easy to solve the form factor equations in this expanded form. For the case rj = —i 
(s = -\) we get 



two 



u(0, 



w 3 (0 = -W-(0, 



(D.56) 



where 



u(0 
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whereas for the case rj = i (s = j) we get 

/-V2\ 

U m (0= 1/2 \U(0, w 3 (H) = -W-(0 (D.58) 




and 



u(0 



ic 



2 cosh | 
3ic (£ — 2iri) 



2 

3ic (£ + 7ri) 
16 cosh | 



As discussed in Section 4, there is a logarithmic piece in the subleading term 
of the 0(4) form factors. As we will see later, such terms are also present in the 
asymptotic expansion of the spin ±1/4 form factors. Anticipating this fact we extend 
(|D.53j) by a logarithmic piece of the form 

^1 ln(A - „) k(* - ffl + + - } . (D.60) 



(A-a)P y ' \ " tv ' ' A- a 

The functions U m and W m satisfy very similar equations to the ones discussed above 
for U m and W m . The main difference is that the residue equations are free of logs. We 
find that non-vanishing solutions are only possible for p = 1/2 or p = 3/2. Finally 
requiring regularity (also at infinity) eliminates all but one possibility: p = 3/2 for 
s = — 1/4 with solution 

U { _-\0=^e-^, (D.61) 

with arbitrary constant £q. This means that the full asymptotic expansion (the sum 
of (ITX531) and (fUBO]) ) is correctly given by (ITX531) alone, with solution (|D.56IID.59j) . 
provided only that we allow the "constant" go depend (linearly) on In A. The actual 
value of go can be calculated from the large A expansion of the exact solution. We 
now turn to this calculation. 

D.4 Contour integral solution 

H. Babujian et al. [18] found the solution of the reduced form factor equations in 
terms of a contour integral: 

3 

fm(Oi,02,O 3 ) = -^ [ due~ 2su t m {e u e 2 ,e 3 ; u) f[<S>(6 , - u) , (D.62) 
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where 



t 1 (9 1 ,e 2 ,e 3 -, u ) 
t 2 {ei,e 2 ,e 3 ;u) 
t 3 {9 x ,9 2 ,9 3 ;u) 
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(D.63) 
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and the contour C (for real rapidities 9i) comes from —00 along a line parallel to 
the real axis and going somewhat below the singular points 6i — iir, then turns back 
and goes around the points 9i before it goes to +00 again parallel to the real line. 
Precisely this integral along such a contour is the special function known as Meijer's 
G-function [19]: 



U 33 



fm(9l,0 2 , O'c 

The parameters depend on the rapidities: 
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Finally we note that Meijer's G-function G33 can be expressed in terms of Gamma 
functions and hyper geometric functions as follows |19j : 



^33 



a\ a-2 a 3 
h b 2 63 



z bl Q(l — a± + bx, 1 — a 2 + 61, 1 — 03 + by, b 2 — bi, 63 — by, —z) + 2 perms . 



(D.66) 



where 



Sl(u 1 ,u 2 ,u 3 ;v 1 ,v 2 ; z) = T (ui)T (u 2 )T (u 3 )T (vi)T (v 2 ) 3 F 2 (u x , u 2 , u 3 ; l-vi,l- v 2 ;z) . 

(D.67) 
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This can be used to express the three-particle form factors as follows. 

2-nis -2se 1 n (l 1 »012 _ 1 _ Wis. iOu itha . 4 ™A 
V 2 ' 2 2?r ' 2 2vr ' 2tt ' 2vr ' / 
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We are interested in the asymptotics of these form factors in the limit 6\ — * +oo. 
The exponential part of the form factor asymptotics, which comes entirely from the 
Gamma functions, is 

e -(i+2«)*i (D.71) 

for the first of the three terms for all f m and is 

e^ ei (D.72) 

for the second and third terms. Thus the exponential part of the asymptotics is 
given by (|D.72jl . which comes from the second and third terms in almost all cases, 
except for s = —1/4, in which case also the first terms contribute. 

To calculate the leading asymptotics of our form factors we will need the asymp- 
totic behavior of the generalized hypergeometric functions 3F2 in the case of some 
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of its parameters large. We will use a simple integral representation [20] of this 
function to establish the asymptotic formulae we need in this calculation. 



3-^2(01, 02,03; h, fo;z) 

r(b 2 ) f l 1 h (°- 73 ) 

= — — r / dt^- 1 (l-t) b '-^- 1 2 F 1 (ai,a 2 ;b 1 ;tz), 

T{a 3 )T(b 2 - a 3 ) J 

which is valid for He(b 2 ) > Re( a 3) > in the range \z\ < 1, but can be extended to 
the limit \z\ — ► 1. The Gauss hypergeometric function 2 Fi can in turn be expressed 
as the integral 

2 F 1 (a, b; c; z) = r(b ^_ fe) £ di t b ~ l (1 - t)^ 1 (1 - tz)~ a , (D.74) 

valid for Re(c) > Re(6) > 0. 

Using both integral representations above simultaneously we can show that for 
large (real) A 



F 2 (ai,a 2 + iX,a 3 + iX; (3% + iX, f3 2 + iX; 1) 



3 

_ + & - «2 - «3 - Ol) 1 A ,a a ^J-sgnA ( D ' 75 ) 

~ r(/3 1+ /3 2 -a 2 -a 3 ) ' ' 

valid for Re(/3 2 ) > Re(a 3 ) > and Re(/3i) > Re(a 2 ) > 0. We are sure however 
that this estimate holds in a larger range of parameters. In particular to get the 
contribution of the first terms for fi, f 2 we need an estimate of the lhs of ()D.75|) for 
ai = l/2,a 2 = i /0 +^,a 3 = -^ + ^,/3i = l + ^,/?2 = |(l + p) + ^for the two 
cases p = ±1. We have numerically checked that the estimate (|D,75P is indeed valid 
for these cases in the range | ^23 1 < 47r. It is plausible that it can also be proved 
for arbitrary values of 3ft 9 23 for some range of 9 9 23 by assuming analyticity of the 
formula in this variable. Applying then (|D.75|) to the first terms of (|D.68HD.70P (for 
s = —1/4) we find 

/W (flljftl ^ )w ^e-^5i__ I 1 ] . ( D.7G, 

^12 

Thus the first terms only contribute to W- for s = —1/4. Their contribution is: 




W[~\0 * ' = -ice~^. (D.77) 
For the second and third terms (of /1 and f 2 ) we have to consider, for large 
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positive A = -j 2 -, the following product: 
r(ai)r(a 2 )r(-ai - a 2 )r(a 3 + i\)T ( - — a<3 — iX 



, 2 ' 1 \ (D - 78 > 

3 F 2 I a 1 ,a 2 ,a 3 + ?A; 1 + a\ + a 2 , - + a 3 -MA; 1 J . 

Its asymptotic form can be established using the integral representation (|D.73P to- 
gether with the formula [21j 
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' + (1 - t) [ln(l - t) - (1) - *(2) + (1 + 01) + (1 + o 2 )] 
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valid in the vicinity of t = 1. With the help of this formula we can calculate the 
large A expansion of (|D.78P : 
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X = *(1 + ai) + (1 + o 2 ) - (1) - *(2) + i — . (D.81) 
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Using (|D.80p in (|D.68p and (|D,69p we find the following results. For s = 1/4 the 
c fy 2 12 terms cancel and only the c terms remain. In this case these contribute 

#12 #12 

to the leading terms U m and we find 



U i+) (0 = -^TT, (D.82) 



2 cosh 
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the same result as we found in the previous subsection. For s = —1/4 (|D.80p gives 
contributions both to the leading U m and the subleading W m terms. We find 

U(~\0 = -P^, (D.83) 

TT q — ITT 

w[-\0 = --e-i^-^ (D.84) 

7T ^ ~~ m 
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and 




(D.85) 



1 T /3\ 1 /l\ T / 1 \ 1 , A i7r 
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2 V 2 / 2 V 2 / 2 2vr 2 



Again, the above results are in agreement with the ones obtained in the previous 
subsection solving the asymptotic form factor equations. Finally we get 
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